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Errata 
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I4.     Poxirth  line  from  bottom:   replace  "  surve"  by  "curve" 
12     Middle  of  page:   replace  "they  are  about  Iw-  feet"  by 

"they  are  also  quite  small  until  the  last  week  when 

they  reach  a  maxim\am  of  1^  feet." 
33     Eighth  line  from  top:   replace  "functions"  by  "function" 
39     Second  line  from  bottom:   should  read  "modification" 
51     Eighth  line  from  top:   insert  "to"  after  the  word  "over" 
5I4.     Middle  of  page:   delete  "over  the  banks"  and  replace  by 

"into  the  river  from  tributaries  and  local  drainage 

between  Wheeling  and  Cincinnati" 
59     Ninth  line  from  footnote:   replace  "solvable"  by 

"solved" 
62     Middle  of  page:   delete  "Cpumps)" 

65  Sixth  line  from  bottom:   replace  "outlines  by 
"outlined" 

66  Eighth  line  from  top:   replace  *4f"  by  "If" 

67  Fourth  line  from  bottom:   repl&be  "tributfeli*lli$"  by 
"turbines" 
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NUI-ERICAL  SOLUTION    OP  FLOOD   PREDICTION  kWD 
RIVER   REGUL/iTION   PROBLEMS 

REPORT   III 

Results  of  the  Numerical  Prediction  of  the  19l|.5  and  19ll-8 
Floods  in  the  Ohio  River,  of  the  19)i7  Flood  through  the 
Junction-  of  the  Ohio  and  Mississippi  Rivers,  and  of  the 
Floods  of  1950  and  19^8  through  Kentucky  Reservoir, 

§1.   Introduction  and  Summary 

In  two  previous  reports  (Reports  I  and  II  from  the  New 
York  University  Institute  of  Mathematical  Sciences,  Nos .  200 
(1953)  and  205  (195ij-))  the  basic  general  theory  for  the 
numerical  analysis  of  flood  wave  problems  in  rivers  was 
developed  and  applied  to  simplified  models  of  the  Ohio  River 
and  of  its  junction  vjith  the  Mississippi  River.   It  xjas  found 
that  the  numerical  treatment  of  flood  problems  seemed  feasible 
from  the  standpoint  of  the  amount  of  calculating  time  needed 
for  a  digital  computer  of  the  type  of  the  UNIVAC  and  thus  it 
seemed  likely  that  problems  for  actual  rivers  could  also  be 
successfully  solved. 

The  present  report  has  as  its  purpose  the  application  of 
the  methods  described  in  the  first  two  reports  to  the  flood 
wave  problems  indicated  in  the  title  of  this  report.   It  might 
be  said  at  the  outset  that  the  results  of  these  calculations 
show  the  numerical  method  to  be  feasible  and  practical. 

The  data  for  the  flood  in  the  Ohio  i^rere  taken  for  the 
case  of  the  big  flood  of  19i|5  and  predictions  were  made 
numerically  for  periods  of  two  weeks  or  more  for  the  375  mile 
long  stretch  of  the  Ohio  extending  from  Wheeling,  West  Virginia 
to  Cincinnati,  Ohio.  For  the  I9I4-8  flood  in  the  Ohio,  6  day 
runs  were  made  during  open  river  conditions.  For  the  flood 
through  the  junction  of  the  Ohio  and  Mississippi  Rivers  the 
data  for  the  1914-7  flood  were  used  and  predictions  were  made  in 


all  three  branches  for  distances  of  roughly  1;0  miles  from  the 
junction  along  each  branch  and  for  periods  up  to  l6  days.   In 
Kentucky  Reservoir,  which  extends  from  Kentucky  Dam  near  the 
mouth  of  the  Tennessee  River  l8t|.  miles  upstream  to  Pickvjick 
Dam,  flood  predictions  were  made  for  the  flood  of  1950  for  a 
period  of  21  days,  and  further  calculations  for  the  I9IJ.8  flood 
for  7  l/h   days  were  made.   In  each  case  the  state  of  the  river 
or  river  system  was  taken  from  the  observed  flood  at  a  certain 
time  t  =  0;  for  subsequent  tinies  the  inflov;  from  tributaries 
and  the  local  runoff  in  the  main  river  valley  were  taken  from 
the  actual  records.  The  differential  equations  which 
characterize  the  flow  of  the  river  \-jeve   then  integrated 
numerically  with  the  use  of  the  UNIVAC  digital  computer  in 
order  to  obtain  the  river  stages  at  futiire  times.  The  time 
required  to  perform  the  calculations  for  the  Ohio  River  on  the 
Ul'IIVAC  is  at  present  one  minute  for  a  one  hour  prediction,  and 
less  than  one-half  minute  for  a  one  hour  prediction  in  the 
other  two  cases.   (We  estimate  that  the  same  calculations  on 
the  IBM  70l|.  coxjld  be  done  in  less  than  l/l5  f^f  the  time 
required  by  the  UNIVAC.)  The  flood  stages  determined  in  this 
way  were  then  compared  with  the  actual  records  of  the  flood. 

The  fact  that  such  flood  wave  problems  in  rivers  can  be 
solved  in  this  way  is,  of  course,  a  matter  of  considerable 
practical  importance  from  various  points  of  viex^r.   Once  the 
basic  data  for  a  river,  river  system,  or  reservoir  have  been 
prepared  and  coded  for  a  calculating  machine,  it  becomes 
possible  to  solve  all  sorts  of  problems  quickly  andX 
inexpensively:   for  example,  the  effect  on  a  flood  wave  of 
damming  a  tributary,  the  relative  meruit s  of  various  schemes 
for  serial  operation  of  dams,  or  the  preparation  of  tables 
showing  the  Influence  of  changing  the  operating  conditions  in 
large  reservoirs,  are  all  problems  vrhich  can  be  successfully 
attacked  numerically.  These  observations  refer  to  the  use  of 
the  UNIVAC  computer.  The  situation  would  be  even  more 
favorable  If  newer  and  faster  computing  machines  such  as  the 


IBM  7OI4.  were  to  be  n.sed.   It  is,  of  course,  also  of  Interest 
to  compare  the  method  of  numerical  computation  with  the  method 
of  using  hydraulic  models;  the  authors  give  their  ideas  and 
opinions  on  this  point  in  section  8  of  this  report. 

Once  it  becomes  clear  that  the  Integration  of  the 
differential  equations  characterizing  flows  and  wave  motions 
in  rivers  can  be  done  accurately  x-jlthout  an  inordinate  amount 
of  expensive  calcula  ting  machine  time  ~  and  this  was  already 
indicated  by  what  was  done  in  Report  II  mentioned  above  -  it 
follows  that  the  success  of  the  method  for  actual  rivers  hinges 
on  the  possibility  of  obtaining  accurate  data  from  vjhich  to 
calculate  coefficients  and  Initial  and  boundary  data  for  the 
differential  equations. 

There  are  four  such  coefficients  (cf.  equations  (2,1) 
of  the  following  section):   the  cross  section  area  A  and  the 
breadth  B  of  the  river,  the  resistance  coefficient  G,  and  the 
inflow  q  from  tributaries  plus  the  local  Inflow  from  the  main 
valley.  The  first  two  quantities  are  purely  geometrical  in 
character  and  could  in  principle  be  determined  from 
topographic  surveys  (as  is,  in  fact,  done  when  hydraulic  models 
are  built).  The  resistance  coefficient  G  is  determined 
empirically  from  records  of  past  floods;  records  of  discharge 
and  stage  along  the  river  are  needed  for  this  purpose.  The 
quantities  A,  B,  and  G  are  all  functions  of  the  location  x 
along  the  river,  and  of  the  stage  H.  The  Inflov;  q  is  assumed 
to  be  a  known  function  of  location  x  and  time  t.   In  the 
problems  to  be  treated  here,  which  were  set  up  to  test  the 
method  of  numerical  calculation,  this  quantity  was  obtained 
from  the  records  of  an  actual  flood;  in  a  prediction  problem 
it  v.'ould  be  obtained  from  knovjledge  of  the  gaged  flows  from 
the  larger  tributaries  plus  estimates  of  the  local  rixnoff 
from  rainfall.   The  quantities  A,  B,  and  G  are  all  fixed 
functions  of  distance  along  the  river  and  stage  for  a  given 
river, -river  system,  or  reservoir.  The  inflow  q  vrill,  of 
co-urse,  vary  from  one  flood  to  another  as  will  the  initial 
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and  boundary  data.   In  the  two  earlier  reports  it  wavS  seen 
that  a  knovrledfe  of  the  state  of  the  river  at  some  initial 
time  (usually  taken  to  be  t  =  0)  is  necessary  in  order  to 
determine  the  flow  uniquely.  This  means  that  both  stage  and 
velocity  must  be  known  initially;  as  a  rule,  the  initial 
velocity  is  determined  by  converting  discharge  data  into 
velocity  data  by  using  cross  section  areas,  and  the  discharges 
in  turn  are  fixed  from  actual  measurements  or  from  rating 
curves  (in  most  cases,  the  latter).   Since  only  finite  lengths 
of  any  system  are  in  question,  it  is  necessary  also  to 
prescribe  boundary  conditions  at- the  upstream  and  downstream 
ends  of  the  system;  for  tranquil,  or  subcritical,  flow  (the 
only  case  dealt  v;ith  here,  since  only  rivers  with  relatively 
low  flox-j  velocities  were  in  question)  it.  is  necessary  to 
prescribe  one  condition  at  each  endpoint,  v/hlch  might  be  stage, 
or  velocity,  or  a  relation  bet^^reen  the  two.   In  general,  the 
discharge  (or,  what  comes  to  the  same  thing,  velocity)  at  the 
upstream  end  should  be  prescribed,  since  that  would  seem  the 
natural  condition  in  practice,  while  at  the  downstream  end  of 
an  open  river  it  seems  reasonable  to  use  an  average  rating 
curve  to. furnish  a  relation  between  stage  and  discharge;  in 
any  case,  a  condition  there  must  be  prescribed  which  has  the 
effect  of  simulating  the  influence  of  the  river  below  the 
endpoint.   In  the  upper  OhioR3.ver,  however,  we  have  prescribed 
stage  at  both  upstream  and  doxmstream  endpoints  (since  our 
basic  object  was  really  to  check  the  general  feasibility  of 
the  numerical  method).   In  the  problem  concerning  the  junction 
of  the  Ohio  and  Mississippi  Rivers  we  have  carried  out  the 
solution  first  by  prescribing  stage  at  the  two  upper  endpoints 
(ice,,  upstream  in  the  Ohio  and  upper  Mississippi),  and 
afterwards  by  using  discharge  data  at  the  upper  ends.   In  both 
cases  a  rating  si;rve  at  the  lower  end  in  the  Mississippi  was 
used  to  furnish  a  boundary  condition.  Good  results  were 
obtained  in  both  cases.   (As  i-rill  be  explained  later,  it  is 
also  necessary  to  impose  appropriate  continuity  conditions  at 


the  junction  itself.)   In  Kentucky  Reservoir,  vjhich  is  closed 
at  both  ends  by  dams,  the  boiuidary  conditions  at  both  ends 
were  formulated  in  terms  of  discharge. 

Once  the  coefficients  A,  B,  G,  and  q  are  known,  and 
initial  and  boundary  data  are  given,  the  differential 
equations  have,  for  all  future  times,  uniquely  determined 
solutions  for  the  stage  H  and  velocity  v  as  functions  of 
distance  x  along  the  river  and  time  t. 

As  was  explained  in  the  earlier  reports,  the  differential 
equations  are  integrated  approximately  by  the  method  of  finite 
differences,  which  yields  values  for  the  stage  H  and  velocity 
v-at  a  discrete  set  of  points  forming  a  rectangular  net  in  the 
x,t-plane.   In  the  Ohio  River  a  net  with- intervals  2^\x  =  10 
miles  in  the  x-direction  along  the  river,  and  ^t  =  9  minutes 
In  the  t-dlrection  was  used,  and  nets  having  approximately 
these  spaclngs  were  also  used  In  the  other  cases.  This  means 
that  the  coefficients  and  the  Initial  data  need  be  knovrn  only 
at  the  net  points,  and  these  quantities  should  consequently 
represent  averages  over  10-mile  stretches,  but  their  values 
are  taken  at  S-rt^ile  intervals  since  a  staggered  net  i-jas  used. 
(Cur  preliminary  vjork  on  simplified  models  of  the  Ohio  River, 
as  discussed  in  Report  II,  indicated  that  these  interval  sizes 
are  sufficiently  small  to  yield  accurate  enough  approximations 
to  the  solutions  of  the  differential  equations.   Hovrever,  the 
simplified  models  were  such  that  the  differential  equations 
had  constant  coeff Icientsj  in  the  actual  cases  the  coefficients 
are  variable  -  so  much  so  that  lO-mlle  Intervals  are  just 
barely  small  enough  to  yield  a  reasonable  approximation.   More 
vjlll  be  said  on  this  point  later  on.)  The  determination  of 
these  coeff icierits  as  averages  over  10-mile  intervals  was  a 
laborious",  difficult,  and  crucial  part  of  our  task.   It  Is  of 
such  importance  that  the  method  of  doing  It  will  be  described 


The  data  characterizing  the  geometrical  and  dynamical 
parameters  for  375  miles  of  the  Ohio  River  required  the 
tabulation  of  1100  constants,  for  example. 


in  detail  in  §3  below.   It  should  be  emphasized  that  this  heavy- 
task  requires  close  cooperation  with  engineers  familiar  vrith 
the  data;  we  were  particularly  fortujiate  in  having  the 
cooperation  of  the  engineers  from  the  Ohio  River  Division  whose 
names  are  given  in  the  acknowledgment. 

Although  we  describe  in  §3  the  methods  used  by  us  in 
converting  the  basic  data  for  a  river  into  data  suitable  for 
our  method  of  numerical  calculation,  it  is  nevertheless  of 
interest  to  indicate  here  in  summary  fashion  how  it  was  done. 
Consider  first  the  resistance  coefficient  G,  which  depends 
physically  upon  the  roughness  and  also  upon  the  nature  of  the 
cross  section  of  the  river  bed.   It  must  be  obtained  from 
records  of  past  floods,  and  it  would  be  very  convenient  for 
this  purpose  to  have  simultaneous  records  of  flood  stages  and 
discharges  at  points  closely  spaced  along  the  river. 
Unf ortiinately,  meaf^urements  of  discharge  are  as  a  rule  available 
only  at  the  ends  of  rather  wide  intervals,  called  reaches  ^   of 
the  order  of  6O-96  miles  in  length,  even  in  the  Ohio  River, 
for  which  the  data  are  more  extensive  than  for  most  rivers  in 
the  United  States.  Thus  an  initial  estimate  for  the  resistance 
coefficient  is  obtained  as  an  average  over  a  distance 
considerably  greater  even  than  the  Interval  size  of  10  miles 
upon  which  our  finite  difference  scheme  is  based,;  a  linear 
interpolation  from  the  midpoints  of  successive  reaches  was  used 
to  fix  the  values  of  G  at  intermediate  net  points.  As  was 
mentioned  above,  the  cross  section  area  A  and  the  breadth  B  of 
the  river  are  geometric  quantities  which  could  be  obtained 
from  topographic  maps.  However,  such  a  procedure  is  extremely 
laborious  and  time  consuming,  and  since  another  equally 
important  empirical  element,  the  resistance  coefficient.  Is 
known  only  as  an  average  over  each  of  the  reaches,  it  seemed 

Ue  had  anticipated  difficulties  because  of  this  fact,  but, 
fortunately  they  did  not  materialize  in  all  of  the  cases.   For 
example,  in  the  junction  problem  and  In  Kentucky  Reservoir, 
the  first  estimates  from  the  basic  data  for  this  coefficient 
were  changed  very  little  subsequently.   However,  quite 
considerable  changes  from  the  first  estimates  were  necessary 
in  the  Ohio  River. 


reasonable  to  make  use  of  an  average  cross  section  area  over 
each  reach  also.   In  fact,  one  of  the  important  aspects  of  the 
results  to  be  reported  here  is  that  it  is  indeed  possible  to 
make  accurate  flood  predictions  by  using  average  cross  section 
areas  in  an  appropriate  way.  Roughly  speaking,  this  was  done 
in  the  Ohio  River  by  analyzing  data  from  past  floods  in  such 
a  way  as  to  obtain  the  storage  volume  in  each  reach  as  a 
function  of  the  stage,  from  which  an  average  cross  section 
area  is  at  once  determined.  At  intermediate  net  points  the 
area  A  was  fixed  by  linear  interpolations  from  the  midpoint 
of  a  reach  to  the  midpoint  of  the  adjacent  one.   In  Kentiicky 
Reservoir,  however,  we  were  supplied  directly  with  storage 
volumes  (obtained  from  topographic  maps)  for  intervals  of 

about  10  miles  in  length.   The  breadth  B  is,  in  principle,  the 

dA 
derivative  rs  of  the  cross  section  area  with  respect  to  stage, 
an. 

It  was  to  be  anticipated  that  this  quantity  would  be  somewhat 
sensitive,  and  this  proved  to  be  the  case;  how  reasonable 
average  values  for  it  vrere  computed  from  the  data  is  perhaps 
best  left  to  the  detailed  description  in  §3.  Thus  the 
quantities  A,  B,  and  G  are  computed  as  numerically  tabulated 
functions  of  stage  at  each  of  the  net  points  along  the  river. 
However,  in  order  to  save  number  storage  capacity  in  the  UNIVAC, 
these .quantities  were  fitted  to  empirical  curves  (quadratic, 
cubic  J  and  hyperbolic  curves  in  different  cases)  with  a  few 
parameters:   the  details  are  of  some  importance,  and  they  also 
are  discussed  in  63.  The  quantity  q(x,t)  which  yields  the 
inflow  data,  is  of  course  taken  directly  from  the  records; 
any  gaged  flows  from  tributaries  were  put  in  at  the  nearest 
interval,  while  the  ungaged  local  drainage  for  a  given  reach 
was  distributed  uniformly  over  the  intervals  in  that  reach. 
It  would  naturally  be  too  much  to  expect  that  a  given 
flood  would  be  accurately  reproduced  by  numerical  integration 
on  the  first  trial.   In  order  to  improve  the  numerical  results 
it  is  in  general  necessary  to  adjust  the  resistance 
coefficient  and  the  average  cross  section  areas  and  breadths. 
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which  are  obtained  initially  as  averages  in  an  appropriate 
sense  over  those  past  floods  for  vrhich  data  are  available,  to 
bring  this  about.   This  is,  in  fact,  what  is  also  done  vfhen 
making  model  studies,  in  which,  however,  all  adjustments  are 
made  through  varying  the  roughness  of  the  model.   In  effect, 
the  observed  flood  is  used  as  a  means  to  correct  first 
estimates  of  the  physical  parameters.   It  had  been  anticipated 
that  very  extensive  changes  in  roughness  coefficients  vjould  be 
necessary  (as  it  is  with  model  studies)  in  order  to  reproduce 
a  given  flood  accurately.  Actually,  it  turned  out  that  the 
first  estimates  of  the  roughness  coefficient  were  quite  good 
in  two  of  the  three  cases,  but  that  the  average  cross  section 
areas  and  breadths  were  in  need  of  revisions,  particularly  in 
some  reaches  of  the  Ohio  River.   In  the  junction  problem,  and 
in  Kentucky  Reservoir,  no  really  extensive  revisions  of  the 
initial  estimates  for  either  the  resistance  or  the  area  and 

■if 

breadth  coefficients  were  necessary  ,  while  both  of  these 
qtiantities  had  to  be  varied  considerably  in  order  to  reproduce 
the  191-1.^  flood  in  the  Ohio  River  with  reasonable  accuracy. 
What  this  means  is  that  the  type  of  basic  data  available  in 
two  of  the  three  cases  treated  by  us  sufficed  to  fix  accurately 
the  geometrical  and  dynamical  parameters  which  govern  the 
flows,  and  that  consequently  the  differential  equations  in 
these  cases  are  adequate  and  correct  formulations  of  the  basic 
laws  which  determine  the  flows  uniquely.   In  the  case  of  the 
Ohio  River,  it  would  be  necessary  to  make  still  more  revisions 
in  the  coefficients  then  we  have  had  opportunity  to  make  so 
far  before  it  would  be  sure  that  the  differential  equations 
mirror  accurately  the  characteristics  of  the  river.   In 


In  Kentucky  Reservoir,  however,  the  mesh  width  of  10  miles 
was  too  large  for  the  finite  difference  scheme  that  worked 
vTell  in  the  other  two  cases.  VJe  therefore  used  a  more  accurate 
finite  difference  scheme  in  the  calculations  for  Kentucky 
Reservoir.   In  fact,  10  miles  is,  we  now  feel,  just  on  the 
borderline  of  what  is  reasonable  as  an  interval  size  for  the 
simple  and  straightforward  ways  in  which  we  approximated 
derivatives  by  difference  quotients. 


particular,  the  use  of  average  cross  section  areas  obtained 
from  past  flood  records  by  balancing  flows  to  obtain  storage 
volumes  is  probably  not  accurate  enough,  and  ought  to  be 
replaced  by  averages  from  topographic  maps. 

In  §[|.  the  results  of  flood  predictions  for  the  19l-!.5  and 
I9I4.8  floods  in  the  375  mile  stretch  of  the  Ohio  River  between 
Wheeling  and  Cincinnati  are  described  in  detail.   In  Pig.  1.1 
we  give  a  graph  showing  a  typical  re suit j  the  graph  shows  the 
observed  and  the  calculated  stages  for  the  19^p  flood  at 
Maysville  for  I3  l/2  days.  The  time  required  on  the  UNIVAC  to 
make  the  computations  was  6  3/1^  hours.  As  one  sees,  the 
agreement  is  generally  good.  The  error  at  the  crest  of  the 
flood  vjas  0.1  foot,  and  the  maximum  error  (late  on  March  6) 
was  1,6  feet.  Upon  going  back  to  the  basic  data,  and  looking 
at  a  map  of  the  drainage  areas,  it  was  observed  that  the 
vingaged  inflow  in  this  reach  which  was  quite  high  for  a  short 
time  on  March  6,  was  mainly  introduced  not  far  above  Maysville; 
on  the  other  hand,  in  making  our  calculations  (as  vie  have 
mentioned  above)  the  ungaged  inflow  was  distributed  over  the 
whole  reach,  and  it  is  thus  not  surprising  that  the  calculated 
hydrograph  is  smoother  than  the  one  actually  observed.   Later 
on,  the  curves  came  together  again. 

Our  process  of  calculation  in  all  three  of  the  problems 
studied  involves,  at  bottom,  the  replacement  of  the  actual 
river  in  all  of  its  complexity  by  a  model  In  which  average 
properties  (average  areas,  resistance  factors,  etc.)  over 
distances  of  varying  length  come  into  play.   It  turns  out  in 
all  three  of  our  cases  that  the  calculated  stages  given  by 
such  models  agree  x/ell  v/ith  the  actual  stages  all  along  the 
river,  but  that  the  discharges  obtained  by  taking  the  product 
of  the  calculated  velocity  and  the  cross  section  area  of  the 
model  at  a  given  point  may  disagree  widely  vjlth  the  local 
discharge  as  observed  in  the  actual  river  -  at  gaging  stations 
at  the  ends  of  reaches,  especially.  The  reason  for  this 
discrepancy  is  that  the  gaging  stations  are  invariably  placed 


::i  •'.-. 


5l8t  H 


P.  9a 


516- 


514- 


512- 


510- 


508- 


506 


504 


502 


500- 


498 


stage 

in 
feet 


STAGE    AT    MAYSVILLE 
DURING     1945    FLOOD 


10 


at  narrow  portions  of  the  river,  and  thus  the  local  areas  and 
resistance  coefficients  at  such  points  differ  widely  from  the 
averages  used  by  us.   It  is  nevertheless  possible  to  obtain 
correct  local  discharges  by  making  an  easy  supplementary 
calculation  which  has  the  effect  of  passing  back  to  the  actual 
river  from  the  model.   The  method  of  doing  so  is  explained  in 
§L\.   in  connection  with  the  problem  of  the  upper  Ohio  River; 
the  same  method  applies  in  the  other  cases,  but  will  not  be 
repeated  in  the  discussion  of  these  cases. 

It  has  already  been  stated  above  that  quite  extensive 
changes  had  to  be  made  in  the  initial  estimates  for  all  of  the 
coefficients  -  resistance  coefficient  G,  area  A,  and  breadth  B  - 
in  the  case  of  the  Ohio  River  in  order  to  reproduce,  as  was 
done  successfully,  the  observed  stages  of  the  19'-|-5  flood.  The 
fact  that  this  was  necessary  is  already  a  strong  indication 
that  the  basic  data  were  inadequate  (or,  perhaps,  not  used  in 
the  best  way  by  us)  to  characterize  the  Ohio  River  accurately  . 
In  fact,  when  the  differential  equations  were  used  sub- 
sequently to  predict  the  stages  in  the  191-1-8  flood,  the  results 
were  not  accui-ate  at  some  (though  not  all)  of  the  gaging 
stations,   (This  x^as  particularly  true  in  the  region  about 
Huntington.   It  has  perhaps  some  significance  to  add  that  we 
are  told  by  the  engineers  of  the  Ohio  River  Division  that  they 
also  have  difficulty  in  getting  their  calcxila  tions  to  check 
in  this  vicinity.)   In  §I|.  the  details  of  the  results  for  the 
I9I4.8  flood  verification,  and  a  discussion  of  the  possible 
causes  for  errors  and  of  possible  ways  to  overcome  them,  is 
given.   In  a  way,  it  was  unfortunate  that  we  began  our  v/ork 
with  the  Ohio  River,  since  we  are  now  convinced  that  the 
problem  presented  by  the  Ohio  River  is  by  far  the  most  difficult 
of  the  three  treated  by  us,  (for  reasons  which  will  be  given  in 
later  sections),  and  that  it  should  be  studied  more  carefully 


There  are  discrepancies  in  the  storage  volumes  of  as  much  as 
20  /6  when  determined  by  different  methods,  and  one  could  have 
legitimate  doubts  with  respect  to  the  accuracy  of  the  inflows 
between  Wheeling  and  Cincinnati,  which  are  a  very  large  part 
of  the  total  flow. 
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than  our  resources  in  time  and  funds  permitted:   with  our 
present  experience  we  feel  that  we  could  attack  the  problem 
in  better  v/ays. 

In  §5  the  calculations  for  the  19'4-7  flood  through  the 
junction  of  the  Ohio  and  Mississippi  Rivers  are  described  and 
analyzed.  Vie   reproduce  here  in  Pigs.  l,2(a)  and  (b)  graphs 
showing  observed  and  computed  stages  at  Cairo,  the  junction, 
and  at  Hickman  In  the  Mississippi  River  below  the  junction. 
Approximately  I4.O  miles  in  each  of  the  three  branches  is 
involved,  and  discharge  was  prescribed  at  the  upper  ends, 
while  a  rating  curve  relating  stage  and  discharge  was  used  as 
a  boundary  condition  at  Hickman,   The  computation  time  on  the 
UNIVAC  for  the  flow  over  a  period  of  20  days  vjas  about  3  1/2 
hours.  As  one  sees,  the  observed  and  calculated  stages  are 
in  very  good  agreement,  with  a  maximum  error  of  about  6  inches 
at  Cairo  and  about  1  foot  at  Hickman.   It  is  to  be  seen  that 
there  is  a  uniform  bias  at  both  Hickman  and  Cairo  in  the  sense 
that  the  observed  stages  are  lower  on  the  rising  part  of  the 
flood  and  higher  on  the  falling  part  than  the  calculated 
stages.  This  is  doubtlessly  the  result  of  using  a  (single- 
valued)  simple  rating  curve  at  Hickman  as  a  basis  for  fixing 
the  relation  between  stage  and  discharge  that  was  used  as  a 
boundary  condition.  Rating  curves  which  depend  on  the  water 
surface  slope  as  parameter  should  perhaps  have  been  used, 
since  that  is  what  is  actually  observed.   Had  this  been  done, 
the  correction  vjould  have  been  in  such  a  direction  as  to 
decrease  the  discrepancy  between  observed  and  calculated 
stages  since  the  actual  rating  curve  relation  would,  for  rising 
stages,  furnish  a  lower  stage  for  a  given  discharge  than  that 
used  in  the  calculations,  and  just  the  reverse  for  falling 
stages.   One  observes  that  the  bias  at  Cairo  is  still  notice- 
able, though  less  in  value  than  at  Hickman.   Only  minor  changes 
were  made  in  the  cross  section  areas  and  resistance  coefficients 
from  the  values  computed  initially  from  the  basic  data. 
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In  §6  the  calculations  for  the  I9I+8  and  1?50  floods 
through  Kentucky  Reservoir  are  described  and  analyzed.   As  uas 
mentioned  above,  the  reservoir  is  13I4.  miles  long  and  extends 
from  Pickwick  Dam  at  the  upstream  end  to  Kentucky  Dam  down- 
stream. Discharge  data  were  used  to  obtain  boundary  conditions 
at  both  ends  of  the  reservoir.   In  Pig.  I.3  results  of  the 
calculations  for  the  1950  flood  are  shox'Tn,   Stages  were 
calculated  for  a  21-day  period  beginning  Jan,  I4.,  and  a  little 
less  than  l\.   hours  of  UNIVAC  time  was  needed  for  the  computation. 
Stages  at  Pickwick  and  Kentucky  Dams"  (the  ends  of  the 
reservoir),  and  at  Perryville  (close  to  the  middle  of  the 
reserx'oir)  are  shown.  The  agreement  between  observed  and 
calculated  stages  is  seen  to  be  very  good.  At  Pickwick  Dam 
and  Perryville  the  errors  are  of  the  oi-der  of  a  few  inches, 
while  at  Kentucky  Reservoir  they  are  about  1  l/2  feet.   One 
sees,  moreover,  that  even  the  minor  variations  in  the  observed 
stages  are  reproduced  faithfully  by  the  calculated  stages. 
These  results  again  were  obtained  without  extensive  changes  in 
resistance,  area,  and  breadth  coefficients  after  their 
determination  from  the  original  data.  A  revision  and  refinement 
of  the  finite  difference  scheme  v:as,  hov;ever,  necessary  (as  was 
mentioned  above),  because  the  mesh  width  of  10  miles  was  too 
large  in  comparison  with  the  rapid  variations  in  cross  section 
area  and  width  vjith  location  along  the  reservoir.   The 
coefficients  used  for  the  1950  flood  were  then  used  in  the 
differential  equations  to  calculate  the  progress  of  the  19i|8 
flood  through  the  reservoir.   In  this  case  the  flood  is 
reproduced  quite  accurately.  The  progress  of  the  19k-Q   flood 
was  calculated  for  7  l/^-  days  only,  since  the  stages  were  then 
higher  than  those  of  the  1950  flood  (hence  it  has  no  meaning 
to  speak  then  of  a  verification),  and  also  the  increase  in 
stage  was  so  rapid  in  some  places  as  to  make  our  finite 


l^ince  discharges  only  are  prescribed  at  these  points,  the 
stages  are  determined  as  part  of  the  solution  of  the 
differential  equations. 
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difference  approximations  too  crude:   either  a  finer  net  should 
be  used,  or  a  niore  refined  scheme  of  approximation  should  bo 
devised. 

Summing  up,  vje  observe  that  the  numerical  method  proved 
to  be  very  successful  for  the  junction  problem  and  for  Kentucky 
Reservoir,  and  partially  so  for  the  Ohio  River.   The  character 
of  the  results  in  the  first  two  cases  -  in  particular,  the  fact 
that  even  minor  variations  in  the  observed  stages  were 
reproduced  with  only  minor  changes  and  adjustments  in  the 
coefficients  -  convinces  us  of  the  feasibility  and  practicality 
of  the  numerical  method.   The  authors  feel  sure  that  the  Ohio 
River  problem  can  also  be  solved  by  numerical  methods  as 
accurately  as  the  basic  data  permits. 

In  §7  ^'0  set  down  some  suggestions,  based  on  our 
experience,  for  modifying  the  numerical  methods.   In  §8  v;e 
discuss  the  relation  between  the  numerical  methods  used  by  us 
for  flood  and  river  regulation  problems  and  the  method  of 
studying  floods  by  means  of  models  of  rivers.  The  viex-js  and 
interpretations  presented  there  are  those  of  the  authors,  and 
do  not  reflect  necessarily  the  views  of  any  others  who  have 
been  associated  with  us  in  this  enterprise.   In  §9  we  describe 
a  method  by  which  rainfall  data  can  be  converted  into  runoff 
data  to  obtain  the  ungaged  local  inflows  in  a  form  suitable  for 
machine  computation;  this  was  carried  out  for  the  Ohio  River. 

Finally,  the  Instructions  to  the  UinVAC  needed  to  solve 
any  of  the  problems  discussed  here  will  be  made  available, 
though  in  a  resti'icted  number  of  copies  because  of  their  bulk. 


Ik 


§2,   Outline  of  the  Numerical  Methods  Used  for  Solving 

Plow  Probleras  in  Rivers 
The  basic  equations  governing  the  flow  in  a  river  are 


BH.  +  (Av)^  =  q 


(2,1) 


^t  ■*■  ^"^x  ■*■  S^^x  =  -Gv|v|  -  2.  V   . 


In  these  equations  H  denotes  the  elevation  of  the  water 
surface  above  sea  level,  v  the  velocity  of  the  flow  in  the 
river,  G  the  resistance  coefficient,  A  the  cross  section 
area,  B  the  vjidth  of  the  river  at  the  vrater  surface,  q  the 
volume  of  inflow  over  the  river  banks  and  from  tributaries 
per  unit  length  and  time.   These  equations  correspond  to 
equations  (2o3)  and  (2»11)  of  Report  I,   However,  the  term 
A.  in  (2o3)  of  the  earlier  report  is  now  vrrjtten  in  the  ob- 
viously equivalent  form  BII.  ,  while  the  coefficient  of  the 
resistance  term  in  (2.11)  of  Report  I  has  been  simplified 
by  writing  it  as  a  function  G(x,H)  vjhlch  is  to  be  empirically 
determined.   The  derivation  of  these  equations,  vrhich  are  in 
any  case  i-jell   l-cnovm,  has  been  given  in  Report  I,   The  nota- 
tion here  is  slightly  different  from  that  of  Report  I;   in 
particular,  it  might  be  noted  that  the  slope  S  of  the  river 
bottom  does  not  appear  explicitly  in  the  second  of  ecji  ations 
(2.1),   S  is  contained  in  the  term  Involving  H   since  II  - 

X  X 

y  -  S,  where  y  is  the  depth  of  the  I'lver.   (It  should  be 
noted  that  H  is  considered  positive  upward,  but  the  slope 
S  is  taken  to  be  positive  although  the  downstream  direction 
is  taken  as  the  positive  x-direction. ) 

The  solution  of  any  concrete  flood  wave  problem  requires 
the  determination  of  H  and  v  as  functions  of  the  location  x 
along  the  river  and  of  the  time  t.   As  V7as  explained  in 
Reports  I  and  II  these  quantities  can  be  obtained  by  integrating 
the  equations  (2,1)  (once  the  coefficients  A,B,G  are  known) 
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provided  that  the  initial  values  for  H  and  v  are  1-cnown  at 
sortie  tine,  say  at  t  =  0,  and  if  in  addition  the  inflou  over 
the  banks  aid  from  tributaries  is  given  as  a  function  of  time: 
that  is,  the  quantity  q  in  the  first  of  equations  (2,1)  should 
be  supposed  Imown  as  a  function  of  x   and  t.   It  is  of  course 
possible  to  deal  only  with  finite  lengths  of  a  given  river 
or  river  system  and  as  a  consequence  boundary  conditions  are 
needed  at  the  end  points.   In  the  problem  of  predicting  a 
flood  through  the  junction  of  the  Ohio  and  Mississippi 
Rivers,  for  example,  boundary  data  were  applied  as  follows: 
The  discharge  was  assuiTied  kixown  above  the  junction  at 
Metropolis  in  the  Ohio  and  at  Thebes  in  the  Upper  Mississi- 
ppi. At  the  dovmstream  end  of  the  Mississippi  at  Hickmai  it 
was  assumed  that  the  relation  betv;een  stage  and  discharge 
was  kno^^In.   In  addition,  at  Cairo,  the  junction  of  the  three 
branches,  it  x-ias  necessary  to  fulfill  transition  conditions: 
these  took  the  form  of  requirements  that  the  stage  in  all 
three  branches  was  the  same  and  that  the  inflows  from  the 
upper  branches  just  balanced  the  outflot^f  into  the  Lower 
Mississippi. 

As  was  described  in  the  preceding  reports,  numerical 
integration  of  the  equations  (2el)  with  the  given  initial 
and  boundary  conditions  is  to  be  performed  by  using  the 
method  of  finite  differences.  However,  before  writing  down 
the  requisite  form-ulas  Involving  finite  differences  there 
is  some  point  in  rewriting  equations  (2,l)  in  the  so-called 
characteristic  form  for  two  reasons:   first  of  all  it  is 
important  to  determine  the  slopes  of  the  characteristics 
in  order  to  fix  a  maximum  safe  interval  for  the  tirae  incre- 
ment /\t,  and  secondly,  this  form  of  the  equation  is  appro- 
priate for  use  in  computations  at  boundary  points,  namely 
at  upstream  and  downstream  ends  of  the  river  stretch  to 
be  investigated. 


"»  1- 
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In  order  to  put  the  equations  (1,1)  in  characteristic 
form  the  second  equation  is  multiplied  by  ±/r~  ^^'^^   added 
to  the  first  equation  thus  yielding  the  system 

(Av)y_+  BH^-q  +y^  |v^  +  vv^  +  gll^  +  Gv  |  v  j  +  ^  |  =  0 
(2.2) 

(Av)y_  +  BH^-q  -/^  -K  ■"■   vv^  +  gH^  +  Gv] v  1  +  ^j  =  0  . 

These  can  in  turn  be  rewritten  in  the  following  form 

where  ^  =  ^l"  +  [v  +  c]  — ,  and  c  =j^   ,   The  quantity  c 
represents  the  speed  of  propagation  of  small  disturbances 
(in  the  caao  of  a  rectangular  channel  c  =ygy,  with  y  the 
depth  of  the  stream)  and  -rr  denotes  differentiation  in  one 
of  the  two  characteristic  directions.   That  is 

(2.1|.).   ^  H  =  V  +  c 

are  the  characteristic  directions  in  the  x,t  plane  in  viiich 
the  derivatives  of  H  and  v  are  taken  in  equation  (2o3)j  and 
the  solution  curves  of  these  first  order  ordinary  differen- 
tial eqiaations  are  called  tVie  characteristics  corresponding 
to  a  given  solution  H(x,t),  v(x,t)  of  (2,1). 

The  method  of  finite  differences  is  based  on  the  deter- 
mination of  approximate  solutions  of  the  differential  equa- 
tions in  a  discrete  net  of  points  in  an  x,t-plane.   There 
are  various  procedures  X'/hich  can  be  used  to  determine  such 
approximate  solutions.   In  general,  use  was  made  of  a 
staggered  rectangular  net  as  Indicated  in  figure  2,1." 


'The  reasons  for  choosing  such  a  net  are  discussed  in  the 
earlier  reports. 


— ...    V        p 
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Figure    Sol 
Staggered  net  point  lattice  in  Interior 

The  general  idea  of  the  method  of  finite  differences  is  to 
advance  the  solution  step-wise  in  time  intervals  of  length 
/\t.      Suppose,  foi'  example,  the  values  of  H  and  v  have 
already  been  obtained  in  a  certain  horizontal  row  of  ret 
points,  say  in  a  row  containing  the  points  L  and  R  of 
figure  2,1,   The  method  of  advancing  the  solution  to  the 
next  row  is  then  as  follows:   Consider  the  point  M  midway 
betx-reen  the  points  L  and  R,   The  valiies  of  H  and  v  at  this 
point  are  defined  as  the  follovring  averages 

(2.5)      Hj^  =  |(Hj^  +  H^)  ,   Vj^  =  \{v^  +   Vj.  )  0 

The  derivatives  of  H  and  v  at  M  are  ap  proximated  by  diffe- 
rence quotients  in  an  obvious  way;   and  these  approximations 
to  the  derivatives  at  M  are  then  inserted  for  the  corres- 
ponding derivatives  in  equations  (2,1 )«   The  result  is  a 
pair  of  algebraic  equations  which  can  be  solved  to  yield 
approximate  values  for  v  and  H  at  point  P,   The  results  are 
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.   ^    [G   v2   +   G   v2]    -   Atf    !R.-   +  !l    ]    q^, 


(2.6) 


Hp  =  I  fH,.  H,]  -  t4W  ^''''"  "  '""  ^ '  '^''  ^^^^ 


The  criterion  for  convergence  of  the  finite  difference 
scheme  as  A^  and  />t  tend  to  zero  is  that  P  should  always 
lie  within  a  triangle  formed  by  the  segment  LR  and  the  two 
characteristics  issuing  from  L  and  R  of  slope  v  +  c, 
given  by  equation  (2ol4.)o 

In  order  to  compute  values  of  stage  or  discharge  for 
net  points  at  either  the  upstream  or  downstream  end,  it 
is  necessary  to  have  given  one  physical  condition,  such 
as  stage  or  discharge  or  a  relation  between  stage  and  dis- 
charge such  as  a  rating  curve.   The  physical  condition  is 
then  used  with  one  of  the  differential  equations  in  charac- 
teristic form  to  determine  the  boundary  values  of  H  and  v. 
For  a  boundary  point  on  the  upstream  side  the  appropriate 
cliaracterlstic  equation"  (i.e«  the  second  equation  of 
(r::,2))when  put  in  finite  difference  form  is  used  first 
with  respect  to  a  point  K,  not  in  the  staggered  lattice, 
and  then  for  the  point  N  which  is  in  the  staggered  lattice, 
as  indicated  in  Fig,  2,2  below*   (we  recall  that  the  posi- 
tive X  -  direction  is  taken  downstream). 


"See  Report  I,  p,30  for  an  explanation  of  this  type  of 
procedure. 


K  (■■■-,   o  P 
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^  X 
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Figure   2,2 
Net  point  scheme  at  left  boundary 

That  is,  we  use  the  given  physical  condition  together  with 
equations  (2c7)  and  (2,8)  belo^^r  to  determine  H„  and  v^   in 
terms  of  already  known  quantities: 


(2,7) 


^'71^     ~"27^5E 'Irl 
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^E-^L 
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"  ^L^L 


=  0 


(2,8) 
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^F-^K 


2   2 

^P-^K 
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^ 


=  0. 


By  qpj_  and  q^^  we  mean  the  inflows  over  the  appropriate  seg- 
ments at  the  appropriate  timoo 

At  the  downstream  end  a  similar  pair  of  equations  (2,9) 
and  (2,10)  is  set  up,  as  indicated  in  Reports  I  and  II,  v;ith 
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the  additional  point  K,  as  indicated  in  Figure  2.3 
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Figure        2,3 
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In  the  case  of  the  junction  problem,  the  values  of  H 
and  V  belonging  to  the  net  point  at  the  junction  are  compu- 
ted from  the  conditions  that  the  stages  in  each  of  the  three 
braiches  are  equal  and  that  the  volume  of  water  I'jhlch  flows 
in  from  the  upstream  branches  leaves  through  the  downstream 
branch,  together  with  equation  (2,9)  and  (2,10)  for  each 
of  the  upstream  branches  and  equations  (2,7)  and  (2,8)  for 
the  downstream  branch  (i.e.  the  lower  Mississippi),   The 


•  ' 
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various  conditions  to  be  satisfied  at  the  junction  can  be 
reduced  to  a  single  implicit  equation  for  the  stage  at  the 
junction  in  the  form 

(2.11)  H  =  F(H) 

where  P  is  an  explicitly  knox^m  function  of  the  unknown  junc- 
tion stage.   Equation  (2,11)  is  solved  by  a  process  of 
iteration,  in  which  the  first  guess  at  the  solution  is 
taken  to  be  the  stage  H  at  the  previous  t5.rae.   This  value 
is  inserted  in  P(H)  to  obtain  a  corrected  value  of  H,  etc.; 
this  process  converges  rapidly. 

In  the  Kentucky  Reservoir  problem  it  was  at  once  no- 
ticed that  the  finite  difference  methods  discussed  above, 
which  were  satisfactory  for  the  Upper  Ohio  and  for  the 
junction  of  the  Ohio  and  the  Mississippi,  failed  to  give 
correct  results.   The  reason  for  this  failure  is  the  rapid 
variation  in  the  coefficients  of  the  differential  equation 
iNTith  the  distance.   In  other  words,  the  essential  quantities 
A  and  B,  the  area  and  width  of  the  reservoir,  varied  m,uch 
too  rapidly'"'  in  relation  to  the  mesh  width  2/\x   of  10  miles, 
to  permit  a  good  approximation  to  the  solution  of  the  dii'fe- 
rcntial  equations;   in  fact,  the  numerical  results  ^^^ei'e  so 
i-7ild  as  to  indicate  strong  divergence.   In  addition,  it 
seems  likely  that  the  difficulty  was  aggravated  bj?  our  use 
of  a  staggered  net:   if  the  variations  in  the  coefficients 
happened  to  be  roughly  periodic  with  a  period  of  ten  miles, 
i.e.  viith  a  period  equal  to  the  mesh  width,  it  is  clear 


"It  should  be  said  that  x^re  vjere  in  a  position  to  ]-:now  these 
variations  in  the  quantities  A  and  B  accurately  since  the 
Tenn.  Valley  Authority  had  furnished  us  with  excellent  data 
in  the  form  of  averages  for  these  quantities  over  10-mile 
intervals.   The  changes  of  i^fidth  and  area  are  sometimes  quite 
abrupt  -  they  vary  by  a  factor  of  tv/o  or  three  in  adjacent 
segments . 
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that  the  shifting  of  the  net  back  and  forth  at  each  time 
step  A*  =  9  Tain,  could  easily  result  in  the  building  up 
of  a  systematic  error.   In  fact,  this  does  happen. 

The  obvious  vjay  to  overcome  the  difficulty  would  be 
to  decrease  the  mesh  width  in  the  x-direction  from  ten 
miles  to  not  more  than,  say,  three  miles.   But  since  the 
time  step  ^t  would  have  to  be  decreased  in  the  same 
proportion  the  calculating  time  on  a  digital  computer 
would  be  increased  by  a  factor  of  9  at  least.  This  ijould 
make  the  use  of  the  Univac  somewhat  impractical,  but  it 
would  not  necessarily  matter  if  a  faster  machine  such  as 
the  IBM  70^  were  to  be  used.  We,  however,  were   using 
the  Univac,  and  hence  found  it  necessary  to  devise  a 
different  way  to  overcome  the  difficulty.   Since  the  diffi- 
culty was  felt  to  arise  because  of  the  shifting  back  and 
forth  of  the  net  in  the  x-direction  at  the  time  intervals 
/\t ,    it  was  thought  that  it  might  be  eliminated  by  approxi- 
mating time  derivatives  through  the  use  of  values  at  the 
time  t  -  /^t   in  addition  to  those  at  the  time  t  in  order 
to  advance  them  to  the  time  t  +  /^\t   -  while  in  all  pre- 
vious schemes  values  at  time  t  only  were  used  as  a  basis 
for  advancing  the  solution  to  the  time  t  +  /jju.      This 
more  complicated  method  of  approximating  time  derivatives 
(which  also  complicates  the  coding  for  the  Univac)  turned 
out  to  be  much  more  accurate,  and  it  led  to  a  satisfactory 
solution  of  the  problem. 

The  new  scheme  for  approximating  derivatives,  v/hich 
we  refer  to  as  a  centered  difference  scheme,  is  described 


'We  have  been  led  to  wonder  whether  some  of  the  difficulties 

encountered  in  the  Upper  Ohio  might  not  also  be  partially 

overcome  by  using  this  teclrinique,   liore  will  be  said  on 
this  point  later. 
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on  the  basis  of  Fig.  2oi^-.   The  essential  new  factor  is  that 
the  values  of  v  and  H 


J.  ^ 
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F  i  g  u  r  e   2«ij. 
Centered  Net  Point  Scheme 

are  advanced  to  the  point  P  by  making  use  of  the  icnovjn 
values  of  these  quantities  not  only  at  points  L  and  R  but 
also  at  point  M  as  v/ell.   In  fact,  the  time  derivatives 
are  calculated  using  values  at  P  and  M,  while  x-derivatives 
are  computed  in  the  sarae  way  as  in  the  stagn;©red  scheme 
described  above,  l^^hlch  requires  using  values  at  L  and  P. 
only.   The  result  is  the  follox-;ing  difference  equations 
as  approximations  to  the  differential  equations  (2,11): 

^P"^II  ^  (^R  Jjil  ,     (V\)  ,  ^   2  ,  qLr7M       ^ 


'M 


(2.11)  and 


^VfMl  ^  (^R^R  -  Vl. 
^M  2A.t        27Yx    ~  "  '^LR 


These  equations  can  be  solved  for  Vp  and  Hp  as  follows 
to  yield  the  approximate  solution  at  point  P  determined 
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from  the  known  values  at  earlier  times: 


(2,12) 

F   =11   -2^y_(_RR  L_Ll  _  { 
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§3.   Methods  of  Obtairing  Coefficients  and  Initial 

and  Boundary  Data  for  the  Differential  liquations 
frora  Empirical  Data.  Fixing  of  the  Maximum 
Permissible  Time  Step, 

In  the  differential  equations  (2.1)  we  must  fix  the 
cross  section  area  A,  the  width  of  the  river   3,  end  'che 
resistance  coefficient   G  all  as  functions  of  the  distance 
X  along  the  river  and  of  the  ste^e   H  (actually  only  at 
points  used  in  the  finite  difference  scheme).   These  co- 
efficients are  fixed  once  for  all  for  a  given  river  no 
matter  what  special  problems  are  to  be  solved,   Tiie  remain- 
ing coefficient,  characterizing  the  flow  over  the  banks 
and  from  the  tributaries  as  fixed  by  the  quantity  q, 
will  differ  from  problem  to  problem  depending  on  the 
known  or  estimated  run-off  frora  rainfall  and  the  inflov/s 
from  tributaries.   It  is  supposed  given  in  any  specific 
flood  wave  problem. 

The  method  of  obtainin.t  the  coefficients  frora  the 
basic  data  depends  upon  the  river  itself  and  also  upon 
the  type  of  information  available.   Actually  several  dif- 
ferent methods  have  been  used  by  us:   at  least  tvjo  for  the 
Ohio  River,  another  for  the  junction  of  the  Ohio  and 
Mississippi,  and  still  another  for  the  Kentucky  Reservoir, 

a)  Determination  of  A  and  B  as  functions  of  x  and  H: 

We  begin  by  discussing  the  determination  of  the  cross 
section  area  A,   This  coefficient,  as  vias   indicated  above, 
must  be  known  as  a  function  of   H  at  every  net  point. 
Ideal  for  this  purpose  would  be  actual  cross  section  areas 
frora  topographical  surveys  at  points  close  enough  so  that 
accurate  average  cross  sections  over  ten  mile  intervals 
could  be  obtained,   (As  vie   have  said  repeatedly  our  net 
points  have  been  chosen  ten  miles  apart  in  the  staggered 
scheme  of  net  points).   Of  course  vjhat  is  wanted  is  an 
average  cross  section  area  for  ten  mile  stretches  centered 
at  each  of  the  net  points.   For  the  case  of  the  Kentucky 
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Reservoir  this  information  vras  furnished  directly  since 
storage  volumes  "for  level  pool"  were  provided.   That  is, 
the  volurae   V   of  the  water  in  ten  mile  stretches  was  given 
as  a  fu.nction  of  the  stage  H^   the  average  cross  section 
area  i;as  then  simply  obtained  by  dividing  the  storage  vol- 
ume by  the  length  of  the  segment.   Such  detailed  data  were  ^ 
not  available  for  the  Ohio  River  between  ITieeling  and 
Cincinnati,  and  other  means  of  determining  average  cross 
section  areas  had  to  be  devised.   It  would  have  been  pos- 
sible and  probably  would  have  given  better  results  to  do 
it  from  topographic  maps;   this  is  v;hat  is  done  in  con- 
structing a  model  of  a  river.   However,  such  a  procedure 
is  extremely  laborious  and  time  consuming  and  it  was  thought 
preferable  to  find  out  to  what  extent  such  refinements  are 
necessary,  or  rather  to  find  out  whether  rougher  and  quicker 
methods  of  determining  cross  sections  would  not  be  just  as 
satisfactory, 

V.e  proceed  to  describe  the  methods  ut;ed  by  us  for  de- 
ternmlng  the  cross  section  area  A  for  the  375  mile 
stretch  of  the  Ohio  r;iver  froia  VJheeling  to  Cincinnati. 
There  aro  main  gaging  stations  at  VJheeling,  3t,  liarys,  Poraeroy, 
Huntington,  Haysvllle,  and  Cincinnati,  as  indicated 
schematically  in  Pig,  3.1.   At  these  stations,  X'jhlch  are 
from  sixty  to  ninety  miles  or  so  apart,  the  cross  section 
areas  are  knoxjn  as  a  function  of  sta  ve.   One  mir-ht  think 
It  reasonable  to  take  the  cross  section  areas  at  the  meas- 
uring stations  as  representative  of  the  channel.   Unfor- 
tunately, these  measuring  stations  aro  located  in  general 
at  points  where  the  river  is  relatively  narrow  (vjl;.ich  is 
natural,  of  course,  since  this  facilitates  measurement  of 
the  discharge)  so  t  :at  these  areas  do  not  represent  a  mean 
croT's  section  over  a  reach.   In  fact  a  comparison  ol'  the 
mean  cross  section  areas  for  a  reach  as  obtained  belov;  with 
the  areas  at  gaging  stations  shows  the  latter  to  be  as  little 
as  half  the  average  area  over  a  reach. 
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Fig.   3.1 
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However,  from  analyses  o£   the  records  of  past  floods 
(which  included  discharge  as  i-jell  as  sta^-e  measurements 
at  the  gaging  stations  raentioned  above)  the  storage  vol- 
umes in  each  reach  were  known  as  a  function  of  the  dis- 
charge at  the  lower  end  of  the  reach.   (The  data  were 
available  in  this  form  because  it  is  in  this  form  that  it 
is  used  in  the  conventional  flood  routing  procedure.)   In 
addition  to  the  storage  volume  as  a  function  of  discharge 
at  the  lovjer  end  of  the  reach  there  is  also  available  the 
rating  curve,  that  is,  the  relation  between  dischr-rge  and 
stage  at  a  gaging  station  for  steady  flow  conditions  in 
the  river.   (Since  actual  flows  are  rarely  steady  some  man- 
ipulation of  the  observational  material  is  necessary  in 
order  to  obtain  a  rating  curve  for  steady  flow  conditions ), 
Wext  it  is  necessary  to  know  the  average  slope  of  the  x^jater 
surface  over  a  reach  for  steady  conditions  as  a  function  of 

the  stage"  at  the  lower  end  of  the  reach  which  can  be 

obtained  from  past  flood  records.   It  is  therefore  possible 
to  calculate  the  stage  at  the  middle  of  the  reach  as  a  fun- 
ction of  the  discharge  at  the  lower  end  of  the  reach,  whence 
the  storave  volume  is  known  for  the  reach  as  a  function  of 
stage  at  the  center  of  the  reach.   One  need  only  divide  the 
storage  volume  by  the  length  of  the  reach  to  define  the 
average  crossection  area  as  a  function  of  stage  at  the  mid- 
point of  the  reach.   This  was  done  foi'>  each  of  the  five 
reaches  in  the  stretch  between  Wheeling  and  Cincinnati, 
Piowever,  as  indicated  above,  this  gives  averages  at  points 
seventy  to  el'^.htj   miles  apart  while  a  cross  section  area  is 
needed  for  the  finite  difference  scheme  at  points  only  5 
miles  apart.   It  should  be  remembered  that  the  mesh  vjidth 
is  ten  miles,  but  ^^Je  use  a  staggered  net.   The  areas  at  the 
net  points  were  obtained  by  linear  interpolation  (as  described 
below)  between  the  midpoints  of  adjacent  reaches. 

In  principle  it  would  be  possible  to  code  the  data  for 

the  cross  section  area  as  a  functj.on  of  stage  and  location 

as  a  two-entry  table.   In  practice  it  is  better,  however  to 

.^.  — .~—    _____^_ 

In  Appendix  II  to  Report  II  it  is  seen  that  the  slope  in  a 

steady  flow  in  a  given  channel  is  related  in  a  unique  way  to 

the  stage. 
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approximate  cross  section  areas  at  the  midpoints  of  the 
reaches  by  explicitly  ■  iven  functi'.ns  involvin;3  a  few  par- 
araeters  so  that  the  digital  co;iputer  can  calculate  the  de- 
rlred  values  frora  such  simple  formulas:   in  this  way  the 
storage  capacity  of  the  digital  computer  is  coaserved.   The 
first  method  of  aporoximatinp;  cross  section  areas  tried  by 
us  was  to  use  the  following  quadratic  approximation  formula 

(3.1)  A(x,H)  =  a(x)[H  -  H^(x)]^  +  A^(x) 

in  ^^^hich  a(x),  H  (x),  and  A  (x)  were  first  evaluated  at 
the  midpoints  x.  of  the  reaches  under  study.   This  formula 

for  the  area  at  the  midpoints  of  the  reaches  was  then  ex- 
tended to  the  intermediate  netpoints  simply  by  interpolating 
linearly  betvjeen  the  midpoints  of  two  nei;-hborin[<:  reaches. 
That  is,  the  quantity  a(x)  for  example  is  assumed  given  as 
follox;s:  ^  _^         ^_^ 

(3.2)  s.{x)   =  -^—^  aix)   +  ~-^^     a(x,  ), 

^1  ^'o    °    ^1  ^o     -■■ 

in  v;hich  x  and  x,  are  midpoints  of  two  nelrrhborina:  reaches 
o      1  ^ 

and  x  locates  any  intermediate  net  point.   Similar  formulas 
were  used  to  fix  the  quantj.ties  H  (x)  and  A  (x)  at  inter- 
mediate points.   For  net  points  which  were  not  located  be- 
tween midpoints  of  two  reaches,  as  in  the  vicinity  of  bound- 
ary points,  extrapolation  formulas  of  the  sa  le  sort  were  used 
with  respect  to  the  two  nearest  reaches,   (In  the  case  of  the 
junction  problem  data  from  reaches  outside  the  boundary  points 
were  available  and  thus  in  this  case  the  boundary  points 
played  no  special  role,   A  slight  modification  in  this  process 
of  fixing  a  cross  section  area  vjas  necessary  at  the  junction, 
in  order  to  make  sure  that  the  sum  of  the  cross  section  areas 
of  the  Upper  Mississippi  and  the  Ohio  equaled  the  cross  sec- 
tion area  in  the  Lower  Mississippi, ) 

Once  the  cross  section  area  has  been  fixed  as  a  function 
of  the  stage,  the  width  B  is  determined  from  the  well  known 
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relation  dA  =  BdH  or  by  setting  B  =  |4  .   Figure  3.2  In- 

dicates  the  actual  average  cross  section  area  for  the  reach 
from  St,  ilarys  to  Pomeroy  and  also  shows  hovj  the  area  is 
approximated  by  the  quadratic  formula  (3.1)  (a  third  curve 

given  by  a  hyperbolic  fornula  is  also  drawn sec  discussion 

below).   Figure  3.3  r^hows  the  corresponding  width  B  de- 
termined by  differentiation  of  the  above  area  curves.   As 
one  sees  -  and  this  is  quite  important  frcra  the  point  of  view 
of  the  later  discussion  -  the  empirical  curve  for  the  width 
B  Is  approximated  by  the  straight  line  obtained  by  differ- 
entiation of  the  parabolic  area  curve.   In  the  upper  Ohio 
this  approximation  to  the  x^fidth  B   is  not  accurate  at  the 
higher  stages,  especially  in  the  reach  Pomeroy-Hunt ington: 
the  width  increases  at  a  much  greater  rate.   The  natui'al 
consequence  was  that  the  river  stages  obtained  by  calcul- 
ation at  Poneroy  were  much  too  high,  as  is  shown  by  Fig*  3»k-t 
vjhich  shows  observed  stages  there  compared  v/ith  those  cal- 
culated using  the  parabolic  approximation  formula  for  areas. 

These  observations  indicate  the  need  for  a  more  accurate 
approximation  formula  for  cross  section  areas  in  at  least 
some  parts  of  the  Uoper  Ohio,   Since  the  linear  variation  in 
width  as  a  function  of  stage  furnished  v/ldths  that  appeared 
to  be  too  small  at  higher  stages,  it  was  thought  better  to 
make  use  of  a  hyperbolic  rather  than  a  parabolic  approximation 
for.;aula  for  cross  section  areas.   The  formula  finally  fixed 
upon  for   A  was  the  following: 

(3.3)  A(x,H)  =  —jfl"  +b(x)[H-h(x)]  +  A^(x), 

and  hence  B  =  dA/dH  is  given  by 

(3.1^)  B(x,H)  =  b(x) a(xj__ 


[H-h(x)]2 

Upon  comparison  with  (3.1)  we  note  that  there  are  nov/  four 
parameters  in  (3.3),  i.e.  a(x),  b(x),  h(x),  and  A  (x)  instead 
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of  the  tliree  parameters  in  (3.1).   The  use  of  these  more 
accurate  formulas  (see  Ti-s.  3.2  and  3.3)  led  to  much 
better  results,  as  Is  shovm  In  Pig.  3.^4-,  where  the  observed 
stages  at  Pomeroy  are  covapared  viith  those  calcul->.ted  using 
the  two  different  formulas  for  cross  section  area. 

For  the  problem  of  the  junction  of  the  Chio  and  the 
Mississippi  Rivers  it  was  found  sufficient  to  use  a  parabolic 

interpolation  formula  for  the  Area that  is,  formula  (3.1) 

v;as  used.   Also,  storage  volumes  for  intervals  of  about  20 
miles  in  length  were  given  as  part  of  the  basic  data,  so 
that  accux'ate  average  cross  sections  viere  known  for  intervals 
of  this  length. 

In  Kentucky  Reservoir  also,  the  average  cross  section 
areas  x-jere  simply  obtained  from  storage  volumes  over  10  mile 
intervals.   These  were  obtained  b'-'-  planimeter  from  topo- 
graphic surveys.   However,  quadratic  approximation  formulas 
for  cross  section  areas  were  not  accurate  enough  in  the  upper 
parts  of  the  reservoir;   in  these  portions  the  following 
cubic  approximation  formula  was  used,  vjith  good  results: 

(3.5)  A(x,E)  =  a(x)[n-h(x)]^  +  c(x)[H-h(x)]  +  A^(x),  with 

a(x),  h(x),  c(x)  and  A  (x)  as  parameters. 

It  should  perhaps  be  emphasized  that  all  of  these  approx- 
imation formulas  for  the  cross  section  A  refer,  in  the  f ii  st 
instance,  to  a  point  at  the  center  of  a  given  reach  (in  the 
Ohio  there  vjere  only  five  such  reaches  of  varying  lengths 
and  in  the  other  cases  the,  reaches  were  approximately  ten 
to  twenty  miles  in  length).   Afterwards  linear  interpolations 
between  midpoints  of  successive  reaches  were  used  to  obtain 
the  corresponding  formulas  at  all  netpoints  in  the  manner 
exemplified  by  formula  (3.2)  for  the  parameter  a(x):   that 
is,  the  linear  interpolation  was  carried  out  v/ith  respect  to 
each  coefficient  separately,  and  not  to  tl:e  formula  as  a  whole i 

"In  the  upper  Ohio,  however,  where  the  hyperbolic  interpolation 
formula  was  used,  it  was  necessary  to  modify  this  process  some- 
what because  the  reaches  v;ere  very  long  and  the  parameters  en- 
ter in  a  nonlinear  vray.   In  some  cases  additional  intermediate 
cross  sections  were  inserted  before  linear  interpolation  v/as 
carried  out. 
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Finally,  It  mif[ht  be  added  that  our  experience  in- 
dicates that  the  approximation  forraulas  sriould  bo  checked 
carefully  at  the  lowest  sta:];es,  -where  they  seem  to  be  sen- 
sitive:  otherwise,  it  can  happen  that  a  formula  which  yields 
reasonable  aoproximatlons  for  the  medium  and  high  stares  may, 
yield  absurd  values  (even  negative  areas)  at  very  low  stages, 
b)   Determination  of  the  resistance  coefficient  G(x,H): 

Ive  proceed  to  discuss  methods  of  deteri-:i:' ning  the  resis- 
tance coefficient   G  v\;l;ich,  ol  course,  is  also  needed  at 
all  of  the  net  points  along  the  river  as  a  function  of  the 
stage   H.   One  way  to  deternine  this  function  would  be  to  cal- 
culate it  frorr.  hydraulic  data,  for  exaraple  froiii  the  formula 

2 

(3.6)  G  =  ^—-TTo   • 

(l.i|9)2p>/-^ 

In  this  formula  n   is  Manning's  roughness  coefficient  and  R 
is  the  hydraulic  radius,   A  raore  direct  approach,  however, 
vjas  preferi'ed  and  equation  (3.6)  was  used  only  as  a  check  in 
order  to  see  xvhetrer  reasonable  roughness  coefficients  re- 
sulted from,  the  empirically  deteriamed  v^alues  of   G, 

The  second  differential  equatioii  of  (2.1)  itself  could 
be  used  as  a  I'leans  to  define   G   if  all  other  quantities  in  the 
equation  were  known.   Portunatelj-,  the  terms   v .  ,  vv^_,  and  §-v 

are  in  general  negligible  in  value  compared  with  the  otl^or 
tvjo  at  any  one  instant  of  time"  so  thct   G   can  be  simply 
computed  from,  the  formula 

(3.7)  G  =  -^  . 

V 

This  formula  rightly  says  that  flows  in  large  rivers  are  very 
nearly  steady  flows  in  which  the  velocity  adjusts  itself  in 
such  a  way  that  the  force  of  gravity  down  the  sloping  bed  of 
the  stream  is  very  closely  balanced  by  friction  and  turbulent 
resistance.   From  past  records  of  floods  the  slope  H,  is 

known,  and  the  velocity  v   can  be  computed  from  the  discharge 
records  and  the  mean  cross  section  area  at  the  middle  o'^  the 

"Their  combined  value  is  almost  always  less  than  1  /  of  the 
other  terms. 
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reach  in  the  manner  described  above.   Formula  (3.7)  then 
furnishes  an  average  value  for   G  as  a  function  of  the  stage 
E  for  each  of  the  reaches.   This  was  the  method  used  to 
fix  the  coefficient   G  foi^  the  upper  Ohio  from  recoi'ds  of 
past  floods.   As  v/ith  the  cross  section  area  A,  it  is  con- 
venient to  represent  the  resistance  coefficient   G  b;/  an 
approxiraation  fornula.   For  the  Ohio  Riv.r  the  follov-inf-: 
formula  was  used: 

(3.8)  G(x,H)  =  a(x)[H-h^(x)]  +  fj^-^j  +  %(x), 

vrhere  a(x),  P(x),  h^(x)  and  G^(::)  are  first  determined  at 

the  midpoints  x.  of  the  reaches  so  tha,t  (3^8)  approximates 
the  average  values  obtained  from  (3.7) c  For  other  values 
of  X,  linear  interpolation  formulas  of  the  type  of  (3.2) 
are  applied  to  the  coefficients  a(x),  P(x),  h  (x),  and  G  (x) 

as  in  the  case  of  the  area  coefficients.   In  fip;ure  (3.5) 
we  plot,  as  an  example,  the  average  resistance  coefficient 
G   for  t]:ie  reach  from  i"omeroy  to  Huntinr;;ton.   (The  resist- 
ance curve  is  afterwards  adjusted  on  the  basl..-  of  trial 
calculations  as  explained  in  sectioii  [j.). 

For  the  case  of  the  junction  problem  and  also  for 
Kentucky  Reservoir  ti\e  determination  of  G  vjas  sim.plified 
because  the  basic  data  furnished  by  the  engineers  included 
the  knowledge  of  what  is  called  the  conveyance  factor  K,  which 
is  related  to  our  coefficient   G  by  the  formula 

2 

(3.9)  G  =  i^  .  ^ 

K^    L 

in  which  L   is  the  length  of  the  reach  for  which  the  con- 
veyance factor   K   is  known.   For  the  case  of  the  junction 
and  Kentuck.y  Reservoir  problem,   K  vjas  kno\m  for  reaches 
approxim.ately  ten  miles  in  length.   The  curves  obtained  from 
(3.9)  were  then  approximated  by  a  parabolic  formula.   At  in- 
termediate net  points  linear  interpolation  was  again  used  to 
fix  the  values  of  G.      ...  .. 
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The  best  checic  on  the  correctness  of  our  coefficients 
is  of  course  obtained  by  checkinfi  the  calculated  stashes  and 
discharges  against  the  observed  quantities  in  actual  floods, 
Hovjever,  in  the  case  of  the  resistance  coefficient   G   a 
check  for  peneral  consistericy  is  pos:.ible  usm;:  formula  (3.6), 
In  this  fornula  the  hydraulic  radius  R  is  a  purely  geomet- 
rical quantity  which  is  easily  computed  from  bhe  known  cross 
section  areas  as  functions  of  the  stage  (as  described  above): 
it  differs,  in  fact,  in  all  of  our  cases  very  little  from  the 
mean  depth  A/B.  Hence  the  formula  (3.6)  makes  it  possible  to 
calculate  Ilanning's  coefficient,  the  roughness  factor  n. 
This  quantity  should  vary  betx^een  ,01  and  .15,  but  for  rivers 
of  the  sort  we  are  interested  in  its  value  should  be  somev;here 
near  ,03.   In  all  of  our  three  cases  reasonable  values  for 
n  viere   obtained,  as  follows: 

Ohio  River    .02  <  n  <  .05  (mostly  .025  <  n  <  .035) 
Junction     ,0l5  <  n  <  ,035 
Kentucky  Reservoir   ,01  <  n  <  ,0i{. 
Also,  the  rou2;hness  coefficient  increased  in  general  with 
stage,  as  it  should. 

The  fact  that  Manning's  coefficient  does  not  vary  a  great 
deal  is  a  fact  that  might  be  used  to  study  flows  in  rivers  for 
which  the  data  needed  to  determii.e  our  coefficient   G   are 
insufficient:   a  first  estimate  of  it  could  at  least  be  ob- 
tained, which  might  then  be  improved  gradually  once  more 
information  on  floods  became  available. 

c)   Preparation  of  initial  and  boundary  data: 

It  has  been  stated  a  number  of  times  that  it  is  necessary 
to  know  the  state  of  a  river  at  some  initial  Instant,  taken 
by  us  to  be  at  t  =  o,  and  this  implies  that  the  stage   H 
and  velocity  v  are  known  along  the  river  initially,   V/ith 
the  stage   H  there  is  no  difficulty:   it  is  directly  given. 
The  determination  of  the  initial  velocity  v   in  general  req- 
uires some  calculation  using  the  basic  data. 
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In  the  case  of  the  Ohio  River  and  also  in  the  case  of 
the  junction  problera  it  was  assuned  that  the  initial  state  x-jas 
near  enough  to  a  steady  flow  that  the  second  equation  of  (2,1) 
could  be  used  to  deternine  a  value  for  v   at  each  net  point, 
now  that  values  for   G  have  been  determined  by  the  raethods 
described  above,   (Some  experiments  were  made  with  the  Univac 
by  using  the  observed  initial  velocities,  arid  it  was  found 
that  res  Its  would  be  influenced  only  slightly  since  the  in- 
itial errors  were  smoothed  out  within  an  hour,)   This  means 
that   V   is  fixed  initially  by  the  calculation  from  the  fox-mula 

(3.10)  v^  +  -g^  +  I  V  =  o, 

vjith  q  the  local  run-off  and  tributary  inflow  for  the  par- 
ticular net  point.   Actually,  the  term  qv/A  is  not  important. 

In  Kentucky  Reservoir  discharge  measurements  are  avail- 
able at  the  ends  of  the  reservoir  only,  i.e.  at  Kentucky  Dam 
and  Pickwick  Dam.   In  this  case  the  velocity  at  the  txw  end 
points  v/as  of  course  fixed  by  the  formula  v  =  Q/A  , 
A  linear  interpolation  of  the  discharge  between  the  two  ends 
then  furnished  the  initial  values  of  the  discharge  at  inter- 
mediate net  points,  and  a  divl3lon  bv  the  local  cross  section 
area  then  furnished  the  Initial  velocities  at  the  intermed- 
iate net  pointn. 

In  the  upoer  Ohio,  the  boundary  conditions  assui'ied  by 
us  were  that  the  stu  'e  was  knoi-n  as  a  function  of  time  at 
V'heeling  and  Cincinnati.   At  Cincinnati,  the  down-stream 
end,  it  i"  of  course  somewhat  artificial  to  prescribe  the 
stage---it  should  rather  be  computed  as  one  of  the  Important 
unknowns  as  was  done  in  the  junction  problem  which  vje  describe 
next.   In  the  problem  of  the  junction  of  the  Ohio  and  Mississ- 
ippi Rivers  the  boundary  data  were  applied  in  two  different 
ways',  with  good  success  in  both  cases.   In  the  first  calcul- 
ation the  stages  at  the  upper  ends  of  the  Ohio  and  the 
Mississippi  were  prescribed,  but  in  the  second  calculation 
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these  were  replaced  by  the  discharges,  which  seems  to  be  the 
more  natural  procedui-e  from  the  practical  point  of  vlevj.   In 
both  cases,  however,  neither  stage  nor  discharge  was  prescribed 
at  Hickman  at  the  lov?er  end  of  the  Mississippi,   Instead, 
the  condition  vjhich  is  natural  from  the  practical  point  of 
viex;;  at  the  lovier  end  point  of  a  portion  of  a  rivex'  which  con- 
tinues open  belovj  that  point  was  used  there  vj.'iich  furnishes 
a  relation  between  sta,';'e  and  dischargee"  In  our  calculations, 
an  average  rating  curve  was  used  that  xvas  the  result  of  obser- 
vations of  past  floods.   In  Kentucky  Reservoir,  which  is  closed 
at  both  ends  by  dams,  the  natural  conditions  were  used,  i.e. 
that  the  discharges  were  assumed  knovm  as  functions  of  the 
time:   these  are  the  physical  quantities  v/hich  are  subject 
to  direct  control, 

d)   Determination  of  the  maximum  permissible  time  step  /\ t : 

As  v/as  explained  in  Report  II,  the  maximuiii  time  step  /\  t 

which  can  be  used  to  advance  the  solution  from  time   t   to 
time   t  +  ^t  is  fixed  by  the  inequality 


(3.11)  At 


A 


/^x 


—  v+c 


< 


in  xvhich  /\  x  Is  the  half  mesh  width  (in  our  staggered  scheme 
of  net  points),   v   is  the  velocity,  and   c   is  the  propagat- 
ion speed  of  wavelets  given  by  the  formula 


"A 


(3.12)         c  =  /|^   =  ,^— , 

with  J        the  mean  depth  of  the  river » 

It  turned  out  that  c  varies  between  the  same  limits 
in  all  thr-ee  of  our  cases,  i.e,  between  ?.0  ft, /sec,,  and  30 
ft, /sec,  corresponding  to  average  depths  between  10  ft,  and 
I|.0  ft,,  vjhile  the  maximum  value  of  v  is  of  the  order  of  5 
or  6  ft, /sec,  (It  is  thus  to  be  noted  that  v;e  are  operating 
alvjays  v/ith  flows  with  velocities  far  under  the  critical  vel- 
ocity).  Since  A.x  in  (3.11)  is  five  miles  in  our  cases,  it 

"It  would  perhaps  be  better  to  take  rating  curves  with  the 
slope  H   figuring  as  a  parameter. 
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is  readily  found  that  /_t  should  be  taken  not  greater  than 
about  10  mir.utes.   We  have  taken  it  to  be  9  minutes  to  be 
on  the  safe  side.   This  is  of  course  a  very  sihort  tine  step 
in  coraparison  v;ith  the  tine.'^,  of  the  order  of  weeks,  for  which 
the  flows  have  been  computed,   Kevertheleas ,  it  is  not  pos- 
sible (at  least  not  without  a  radical  revision  in  the  whole 

method  of  computation)  to  relax  this  condition which,  it 

will  be  remembered,  is  imporjed  by  basic  mathematical  fact:;  con- 
cerning our  differential  equations.   In  fact,  if  /_t  is  taken 
even  slightly  larger  than  the  lirait  iraposed  by  (3,11),  the  re- 
sult is  likely  to  be,  not  simply  inaccurate  results,  but  rath- 
er values  of  the  unknowns  which,  oscillate  wildly.   This  ex- 
perience  i-jell  known  in  similar  problems  in  other  fields, 

particularly  in  ^jas  djmamics was  veidfied  for  our  problems 

by  an  empirical  test  usin;;  an  exact  solution  of  the  differen- 
tial equations  as  a  basis  for  comparison.   A  steady  progress- 
ing v.'ave  in  a  uniform  channel  (in  fact,  the  wave  described 
in  Report  II)  was  tal-.en  and  initial  values  were  chosen  to  con- 
form to  it;   these  x^/ere  then  advanced  in  the  time  by  numerical 
calculation  usin~  our  methods,   i/hen  /'.  t  vjas  chosen  properly, 

the  results  checked  the  knovrn  solution  very  v/ell;   hov;ever, 
v;hon  /\ t  was  taken  larrer  than  the  permitted  value  for  conver- 
gence, vjild  oscillations  occurred  at  once  in  the  vicinity  of 
thoce  places  where  condition  (3.11)  vjas  violated.   The  same 
experience  was  noted  also  in  a  simple  model  of  the  Kentucky 
Reservoir, 

e)  Further  remarks  about  the  preparation  of  data: 

In  this  section  x-'e  have  de.'";cribed  hov;  the  basic  data  for 
a  river,  river  system,  or  a  large  reservoir  can  be  used  to 
furnish  the  coefficients  and  other  data  needed  to  formulate  the 
flovj  problems  in  terras  of  differential  equations.   It  should 
be  remarked  that  the  basic  data  for  the  most  part  are  worked 
out  for  steady  conditions,  and  that  such  data  may  not  be  suf- 
ficient in  all  respects  for  the  purposes  in  view.   The  acid 
test  is  to  check  the  results  for  actual  floods  aoialnst  those 
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obtained  by  calculation.   If  the  tv;o  do  not  chock,  this  im- 
pllss  that  there  are  errors  in  cross  s-ction  areas,  or  rej- 
istance  coefficient,  or  inflows,  etc.   These  quantities 
should  then  be  changed  in  such  a  v;ay  as  to  sive  better  re- 
sults.  It  turned  out  in  the  tl-irec  problems  treated  by  us 
that  only  rainor  revisions  were  necessary  in  two  of  the  cases, 
but  that  extensive  revisions  in  a  third  (the  uoper  Ohio  F'.ivcr) 
were  needed.   Actually,  the  method  of  numerical  calculation 
has  in  it  the  inherent  possibility  of  improving  the  bar:ic 
data  for  a  river  by  constant  checking  against  the  rei;ults 
for  ne^^;  floods,  and  making  changes  in  coefficients  where 
changes  are  indicated.   In  the  next  sections,  where  our  fin- 
al results  for  the  ti^ree  problems  are  prer;ented,  some  details 
on  this  matter  \-iill   be  Fiven, 
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il\..      Comparison  of  Predicted  and  Observed  Stages  in  the  Upper 
Ohio  River  in  the  1914-5  and  1914-8  Floods 

In  Fig.  3.1  of  the  preceding  section  a  diagrammatic 
sketch  of  the  Ohio  River  between  Wheeling  and  Cincinnati  is 
shov:n  together  v:ith  the  reaches  and  gaging  stations.   By  the 
methods  discussed  in  the  preceding  section  resistance 
coefficients  and  cross  section  areas  representing  averages 
over  each  reach  are  available.  These  in  turn  furnish  by  linear 
interpolation  values  at  the  net  points  of  our  finite  difference 
scheme.  An  interval  between  net  points  of  10  miles  in  the 
staggered  scheme  described  in  §2  was  taken,  on  the  basis  of 
the  experience  with  a  simplified  model  of  the  Ohio  River  which 
was  presented  in  Report  II.  Calculations  for  the  actual  Ohio 
River  for  a  limited  period  of  time  vjere  made  with  a  5-niile 
interval,  in  order  to  get  some  idea  of  the  possible  errors  to 
be  expected  from  use  of  10-mile  intervals.  The  results  at 
Pomeroy  for  a  36-hour  period  are  shown  in  Fig.  l+.l.   There  is 
a  difference  of  6  inches  betvjeen  the  two,  as  i-ie   see  -  which  is 
not  entirely  negligible.   However,  it  was  nevertheless  decided 
to  proceed  on  the  basis  of  10-mile  intervals  in  order  not  to 
use  too  much  calculating  machine  time."  A  time  interval  of  9 
minutes  was  used  which  is  well  under  the  maximum  permissible 
for  convergence  of  the  finite  difference  scheme. 

Flood  calculations  for  the  191-1-5  flood  were  begun  at  a 
time  when  the  river  was  low.   Calculations  vrere  first  made  for 
a  36-hour  period  during  which  the  flood  was  rising;  as  stated 
earlier  these  calculations  viere  made  using  the  measured  inflows 
of  tributaries  and  the  estimated  runoff  in  the  main  valley. 
Upon  comparison  with  the  actual  records  it  was  found  that  the 


'it  should  be  repeated  here  that  we  began  the  investigation 
which  forms  the  subject  of  this  report  with  the  problem  of  the 
Ohio  River,  and  had  therefore  to  proceed  without  the  aid  of 
previous  experience  from  any  source.   Realizing  that  it  would 
doubtlessly  be  necessary  to  experiment  i-Jith  and  to  revise  the 
methods  of  calculation,  it  vjas  decided  that  to  use  too  fine  a 
net  at  the  beginning  would  result  in  vjaste  of  expensive 
calculating  machine  time. 
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Calculated    Stages     for     Pomeroy 

a-  using    5   mile    intervals 
b-  using    10    mile     intervals 
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calculated  flood  stapes  were  5.n  general  systeraatlcally  higher 
than  the  observed  stages  and  that  the  discrepancies  Increased 
steadily  with  increase  in  the  time. 

For  example,  Fig.  l\.,2a   shows  the  calculated  versus  the 
observed  stages  for  the  first  36  hours  at  Pomeroy.   It  was 
reasonable  to  suppose  that  the  deviation  was  probably  due  to 
an  inaccuracy  in  the  resistance  coefficient.   Consequently  a 
series  of  flow  calculations  was  made  on  the  UNIVAC  in  which 
this  coefficient  was  varied;  from  these  results  an  adjtisted 
resistance  coefficient  was  estimated  for  each  of  the  reaches. 
In  Fig.  [|..2b  the  results  of  two  such  computations  are  shown. 
In  both  cases,  the  resistance  coefficient  G  was  lowered 
substantially  (a  maximum  of  20  %  -30  %  )  over  a  portion  of  the 
river,  follovjed  by  a  recalculation  of  the  flow  for  a  2l4.-hour 
period.  The  curves  show  the  difference  in  stages  obtained 
when  the  resistance  coefficient  is  changed.   In  case  (i),  the 
coefficient  G  was  decreased  considerably  at  and  near  Pomeroy, 
and  the  resialt  was  a  change  of  stage  of  about  0,3  foot.   In 
case  (ii),  the  change  in  G  was  greater,  and  it  extended  over 
a  larger  portion  of  the  river;  the  result  was  a  much  greater 
change  in  stage  -  a  maximum  of  1  l/2  feet  -  as  was  to  be 
expected.   In  both  cases  one  observes  that  the  maximum  lowering 
of  the  stage  occurs  somewhat  upstream  from  the  region  where  G 
is  lowered,  and  is  follov/ed  dovmstream  with  a  somev/hat  smaller 
increase  in  stage.   Lowering  the  resistance  in  one  section 
seems  to  increase  the  floi-j  above  that  section  and  to  pile  up 
the  water  downstream.   Once  the  effects  of  changes  in  G  have 
been  estimated,  it  becomes  poss.ible  to  make  the  changes  in 
such  a  manner  as  to  bring  calculated  stages  into  agreement 
with  observed  stages.  Actually  this  was  done  rather  roughly, 
simply  by  shifting  bodily  the  original  curves  for  the 
coefficient  G  as  a  function  of  stage  by  a  displacement  of  an 
appropriate,  amount  in  each  reach  (i.e.,  only  the  constant 
G  (x)  in  (3»8)  was  changed).   In  other  words  no  attempt  was 
made  to  make  corrections  that  vrould  require  modificatio  in  the 
shape  of  these  curves  in  their  dependence  on  the  stage. 
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THE    CHANGE  IN  STAGE,   AFTER  24  HOURS, 
WHICH  RESULTS  FROM   A  CHANGE    IN 
RESISTANCE     COEFFICIENT 


I^G^-G, 

40x20-« 
ft. 


case  (i) 


A  ^2"^! 


20 
Wheeling 


60  80 

Pomeroy 

Fig.  4.2b 


40         160 
Cincinnati 


ko 


In  Fig.  l4.,3  we  plot  the  average  resistance  for  the  i^each 
from  Pomeroy  to  Huntington.  We  rhow  the  curve  determined  by 
the  engineers  from  the  basic  data,  the  curve  we  foiond  using 
formula  (3.7),  and  the  shifted  curve  determined  by  the  trial 
calculations.   In  Pig.  l4..2a  we  have  already  given  the  stages 
calculated  with  the  original  and  the  shifted  resistance 
coefficients  at  Pomeroy  for  the  first  36  hours.  The  new 
coefficients  thus  corrected  on  the  basis  of  36  hour  predictions 
(and  thus  for  flood  stages  far  under  the  maximum)  were  then 
used  to  continue  calculation  of  the  flow  for  various  6  day 
periods  as  well  as  one  16  day  period,  with  results  to  be 
discussed  in  a  moment. 

It  should  be  said  at  this  point  that  making  such  a 
correction  of  the  resistance  coefficient  on  the  basis  of 
comparison  of  results  from  a  calculation  of  an  actual  flood 
with  the  observed  quantities  corresponds  exactly  to  what  is 
done  in  making  model  studies.   Indeed,  in  making  model  studies 
no  first  estimate  for  the  resistance  is  possible  a  priori  as 
is  the  case  with  the  method  being  described  here;  instead  it 
is  always  necessary  to  make  a  number  of  verification  runs  after 
the  model  is  built  in  order  to  compare  the  flood  stages  given 
by  the  model  with  actual  floods.   In  doing  so  the  first  ri.in  is 
normally  made  without  making  any  effort  to  have  the  resistance 
correct.   In  fact,  the  roughness  of  the  concrete  of  the  model 
furnishes  the  only  resistance  at  the  start.   Of  course,  it  is 
then  observed  that  the  flood  stages  are  too  low  compared  with 
an  actual  flood  because  the  water  runs  off  too  fast.  Brass 
knobs  are  then  screwed  into  the  bed  of  the  model  and  wire 
screen  is  placed  at  other  parts  of  the  model  to  roughen  it 
until  it  is  foiind  that  the  flood  stages  given  by  the  model 
agree  with  the  observations.  This  is  in  effect  what  was  done 
in  making  numerical  calculations  except  that  the  empirical 
data  furnished  at  least  a  first  estimate  for  the  friction 
resistance  in  the  river  channel  (for  the  junction  and  Kentucky 
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Reservoir  problems,  the  data  on  resi.stance  and  area  v/ere  given 
for  shorter  reaches,  and  little  adjustment  of  the  resistance 
coefficients  was  necessary  -  see  §§5  and  6). 

In  Fig.  i^.i^-  the  rer.uDts  of  the  first  calculation  of  the 
I9I4.5  flood  are  indicated.   Shovn  in  the  figure  are  the  river 
stages  at  Pomeroy.  As  v/as  mentioned  above,  the  resistance 
coefficients  used  were  those  obtained  upon  correction  after 
making  36  hour  trial  runs.   The  cross  section  areas  used  v;ere 
those  obtained  by  the  parabolic  interpolation  formula  discussed 
in  the  preceding  section;  and  as  v;as  already  stated  there, 
they  proved  to  be  not  sufficiently  accurate  at  the  hig;hest 
stages.   One  sees  that  the  predicted  flood  stages  approximate 
the  observed  stages  with  errors  at  the  higher  stages  of  some- 
what less  than  a  foot  for  the  first  6  day  period.  At  the 
crest  of  the  flood,  however,  (on  about  March  9th,  that  is) 
after  11  days  the  error  at  the  crest  of  the  flood  is  more  than 
ll-  feet.  As  the  flood  recedes  and  stages  become  less  the  error 
once  more  becomes  relatively  small.   That  the  calculated  flood 
stages  come  out  too  high  at  the  higher  stages  was  judged  to 
be  due,  as  \<!as   explained  above,  to  the  fact  that  the  actual 
cross  section  areas  were  apparently  not  well  approximated  at 
high  stages  bj   the  quadratic  formulas  used  by  us.  We  there- 
fore replaced  the  quadratic  formulas  by  hyperbolic  formulas 
and  obtained  much  better  resuD.ts  as  can  be  seen  in  Pig.  Ii-.i|.. 

Upon  using  the  hyperbolic  interpolation  formulas  for 
cross  section  areas  (again  see  the  preceding  section  for 
details),  the  results  for  the  19U5  flood  were  very  much 
improved  at  all  of  the  stations.   In  general,  the  errors  in 
stage  at  the  crest  of  the  Hood  are  small,  of  the  order  of  a 
foot  or  less,  and  the  maximum  errors  are  of  the  order  of  2   feet 
or  less,  as  one  sees  from  Pigs.  i|.»5a>  b,  c,  d  which  show  the 
river  stages  at  St.  Marys,  Pomeroy,  Huntington,  and  Kaysville. 
In  addition,  the  minor  variations  in  the  observed  stages  are 
reproduced  rather  well  by  the  calculated  stages,  though  the 
latter  in  general  tend  to  be  smoother  curves.  This  smoothing 


583 

581— i 
579- 
577—1 
575- 
573- 

571- 

569- 

567H 


H 

stag' 
in 
feet 


561- 


noon 


STAGE  AT  POMEROY  DURING  1945  FLOOD 


I 


I    I    M    M    '    M    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    M 
27  28  I  2  3  4  5  S 

FEB  MAR 

Fig.  4.4 


I    I    M    '    I    '    I    M    '   I    '    I    M    '    I    I    I    M    '    I    '    I   '    '    '    I       ^ 
7  8  9  10  II  12  13  14 


>^  H 


622- 


stage 

i  n 
feet 


620  — 


618- 


616- 


614 


612  — 


610  — 


608  — 


606 


604  — 


603 


noon  .   noon 


'    I    '    '    '    I    ' 
/      27  28 

(602.85)   FEB 


/^  H 


p.   41  b 


622- 


620  — 


STAGE    AT    ST,  MARYS     DURING 
1945     FLOOD 


618 


616- 


614  — 


612- 


610  — 


608- 


606- 


604- 


603 


/      27 

(602.85)   FEB 


I  2 

MAR 


57  3 


571 


>^    H 
i  n 
feet 


569- 


567- 


565- 


563- 


561- 


559- 


557- 


555- 


553 


noon       no or 

■    I    I    I    I    I 


p.  41c 


57  3 


571- 


569- 


567- 


565- 


563- 


561- 


559 


557- 


555 


553 


STAGE  AT  POMEROY  DURING  1945  FLOOD 


n  I  I  I  I — I — I — I — I — I — I — I — I — 1 — I — I — r 
4      5      6      7 
Fig.  4.5b 


550^^  H 
in 


feet 


548- 


546- 


544- 


542- 


540- 


538- 


536- 


53  4- 


532 


530 


noon  noon 

— I    '    I  •    I   ' 

27  28 

FEB 


550^  H 
in 


feet 


STAGE    AT    HUNTINGTON    DURING 
1945     FLOOD 


P.  4ld 


548- 


546- 


544- 


542 


540- 


536- 


53  4- 


532- 


FEB 


MAR 


Fig.  4.5c 


5I8>^ 

516- 

514- 

512- 

510 

508- 

506- 

stage 
in 

feet 


504- 


502- 


500- 


^98 


p.    41  • 


518 


516- 


514- 


512- 


510- 


508- 


506- 


504 


502- 


500 


^98 


"1    \   1    I   I    I   I   I    I    rn   I    r 
27  28  12 

FEB  MAR 


I    I    I    I     I    I    I    I    I    I    I     I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I 
3  4  5  678  9I0III2 

Fig.  4.5d 


k2. 


effect,  is,  In  sorae  instances  at  least  due  to  the  distribution 
of  the  local  inflow,  which  \ms   done  more  smoothly  in  the 
calculation  than  is  true  in  actuality  (see  the  discussion  in 
the  introduction  with  reference  to  the  hydrograph  at  Maysville, 
particularly  near  March  6). 

By  and  large,  it  is  rip-ht  to  say  that  the  calculations 
reproduced  the  observed  stajres  of  the  19i|-5  flood  with  quite 
reasonable  accuracy  after  a  good  deal  of  adjustment  of  the 
cross  section  areas  and  resistance  coefficients  had  been  made. 
The  necessity  for  such  revisions  of  the  coefficients  obtained 
from  the  original  data  is  a  strong  indication  that  the  original 
data  were  not  accurate  enough  -  especially  near  Huntington. 
That  this  is  true  is  then  borne  out  by  the  results  for  the 
I9I+8  flood  in  the  Ohio.   If  the  coefficients  used  for  the  19l|5 
flood  were  correct  it  should  be  nossible  to  calculate  stages 
correctly  for  any  other  flood.   However,  \^pon  doing  so  for  the 
I9I1.8  flood  the  results  are  not  as  accurate  as  they  should  be. 
Figure  li.Sa,  b,  c,  d  shoxjs  calcula-ued  versus  observed  stages 
for  the  I9IJ.8  flood.  All  of  these  show  the  calculations  for  a 
6  day  period  starting  April  1^ .     As  one  sees,  the  errors  are 
high,  particularly  at  Hujitlngton,  where  the  crest  is  wrong  by 
about  1  1/2  feet.  At  other  stations  the  results  are  better. 
Figure  I;.. 7  shows  the  stages  at  St.  Marys  for  a  longer  period 
starting  April  12  :   here  the  results  are  quite  good.  Much 
more  vjork  would  be  required  to  obtain  more  accurate  coefficients 
for  the  Ohio,  but  there  is  no  doubt  that  they  could  be  obtained. 
In  addition,  various  revisions  in  computational  methods  should 
be  made,  on  the  basis  of  experiences  vrith  the  computations  for 
Kentucky  Reservoir.  Further  discixssion  of  these  points  is 
postponed  to  §7,  after  the  results  for  the  junction  problem 
and  Kentucky  Reservoir  have  been  presented,  and  comparisons 
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The  reason  for  the  shorter  period  at  the  other  stations  is 
that  the  river  was  not  open  throughout  the  period  in  question 
-  instead,  a  dam  at  Gallipolls,  about  120  miles  below  St.  Marys, 
was  in  operation,  and  this  has  considerable  effects.   Since 
St.  Marys  was  rather  far  upstream  it  was  thought  reasonable  to 
ignore  the  effect  of  the  dam  on  the  stages  there. 
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and  contrasts  can  be  made.  Also,  In  §7  suggestions  of  various 
kinds  are  given  for  improving  the  methods  of  computation,  based 
on  our  experience . 

So  far,  our  discussion  of  results  has  centered  around 
the  comparison  of  observed  with  calculated  stages,  v;hile 
velocities  and  discharges  have  received  little  or  no  attention. 
In  general,  the  stage  is  probably  the  quantity  of  most  interest, 
but,  as  we  knov;,  it  cannot  be  computed  by  our  m^ethods  vrlthout 
also  computing  the  velocity  at  the  same  times  and  places.   In 
addition,  the  discharge  is  in  many  cases  the  natural  quantity 
to  prescribe  as  a  boundary  condition.   Our  method  requires  us 
to  replace  the  actual  river  cross  sections  by  averages  -  in 
the  case  of  the  Ohio  River,  by  averages  over  quite  long 
distances  -,  and  the  velocities  we  compute  are  therefore  also 
certain  averages  with  respect  to  distance.   In  all  three  cases 
treated  by  us  the  stages  obtained  from  our  model  of  the  actual 
rivers,  or  the  reservoir,  could  be  taken  as  the  stages  in  the 
actual  river  with  good  accuracy'';  but  if  one  were  to  calculate 
discharges  by  multiplying  the  cross  section  areas  used  by  us 
by  the  velocity  computed  by  us,  the  result  would  often  differ 
very  greatly  from  the  observed  discharge  at  a  given  point. 
In  other  words,  it  is  necessary,  in  calculating  discharges, 
to  make  a  supplementary  calculation  in  order  to  pass  back  from 
our  model  to  the  actual  river.  We  proceed  to  give  a  simple 
way  to  obtain  correct  discharges. 

VJe  know  that  the  velocitj''  at  a  particular  place  at  a 
.given  time  is  in  general  given  quite  accurately  by  the  formula 
(of.  §3): 

(^.1)  Gv^  =  gli^   . 

However,  that  will  be  true  only  if  G,  the  resistance 
coefficient,  has  the  correct  local  value.  VJe  have  used  the 
formula  in  order  to  compute  an  average  value  for  G  over  the 
various  reaches.   How  widely  the  average  value  departs  from 
the  local  valiae  at  the  gaging  stations  can  be  seen  from 


'•/   «•:  ji 


Pig,  L|..8,  which  shows  the  average  values  of  G  for  all  of  the 

reaches  together  ^^;ith  local  values  at  Wheeling  and  Maysville. 

The  local  values  viere  obtained  from  the  local  values  of  the 

slope  H   (these  in  turn  were  knox^rn  from  the  basic  data,  which 

included  stage  measurements  at  two  points  at  the  ends  of  the 

reaches  which  were  only  a  mile  or  so  apart),  and  the  local 

value  of  the  velocity  (from  Q/A,  with  A  the  actual  cross 

section  area).   The  wide  divergence  of  the  average  from  the 

local  values  of  G  is  doubtlessly  the  result  of  the  fact  that 

the  gaging  stations  are  placed  at  exceptionally  constricted 

portions  of  the  river. 

Nevertheless,  correct  values  for  the  actual  local 

discharge  can  be  obtained  from  our  results,  in  the  following 

vjay.   Let  us  introduce  first  a  few  special  notations  for  this 

purpose.  By  A  and  v  we  mean  the  ci'oss  section  area  used  by 
c       c 

us  in  the  computation  and  the  velocity  obtained  from  our 
calculations;  by  A  and  v  \oe   mean  the  actual,  local,  values 
of  the  same  quantities.   This  means  that  the  discharge  Q, 
obtained  directly  from  our  results  XNrould  be  given  by  Q^  -   A^v^, 
vrhile  the  actual  discharge  Q^  is  found  from  Q^  =  A^v^ .  We 
introduce  also  the  symbols  G  and  G^  to  distinguisli  between 
calculated  (i.ec  from  our  averages  over  reaches)  and  actual 
local  values  of  the  resistance  coefficient.   Since  in  all  of 
our  cases  the  stage  H  and  slope  H  are  correctly  given  on  the 
whole  by  our  calculations,  it  follows  from  formula  (l+.l) 
>jritten  down  above  that  the  following  relation  holds: 

()i.2)  G  v'^  =  Gv^      . 

^  *^   '  c  c    a  a 

Since  Q  =  A  v  and  Q„  =  A  v  ,  the  following  relation  for  Q„ 
c    G  c      a    a  a  "  a 

therefore  holds: 

A  fcT 


(J+c3)  \  =  r-jG^^ 

c  <     s. 


c 
c  <  a 


Ue  repeat  the  significance  of  the  terms  in  this  equation: 

A  and  G  are  the  actual  local  values  of  area  and  resistance 
a      a 
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coefficients, to  be  found  from  the  basic  data;  while  A   and  G 

c      c 

are  the  values  used  in  the  machine  computation  for  these 

quantities,  and  Q  =  A  v  with  v   the  value  for  v  obtained  by 

c    c  c       o 

our  computations. 

Figure  Ij..9a,  b  shows  the  results  of  such  a  computation 
at  Wheeling  and  Maysville.  As  one  sees,  the  calculated 
discharges  mirror  the  observed  discharges  very  ^^7ell,  once  the 
discharges  are  obtained  by  using  equation  (!4..3)»  That  the 
results  would  be  quite  wrong  if  Q.  were  to  be  taken  for  the 
discharge  instead  of  Q  is  also  shovm  by  Fig.  )-|..9b,  which 

contains  a  graph  of  the  factor  A„/A^  t^'T/'^'^''      as  one  sees, 
this  factor  applied  to  Q  can  change  the  discharge  by  as  much 
as  100  96  . 

Finally,  we  give  the  result  of  a  calculation  v/hlch 
displays  the  flexibility  of  the  numerical  method.   Once  the 
basic  data  for  the  river  have  been  coded  it  becomes  possible 
to  experiment  in  many  ways  with  respect  to  hypothetical  flows 
which  would  result  from  varying  conditions  in  the  river.  As 
a  case  in  point  we  found  it  very  simple  to  carry  through  a 
computation  in  which  the  flow  from  one  of  the  major  tributaries 
of  the  Ohio  River,  the  Kanawha  River,  was  cut  off  for  36  hours 
-  as  might  be  the  case  if  there  vieve   a  dam  in  the  Kanawha 
River.  The  upstream  and  downstream  effects  of  such  an 
operation  were  easily  evaluated.   In  Fig.  [(..lOa,  b  the  river 
stages  with  and  without  the  flow  from  the  Kanawha  River  are 
shovm  at  the  nearest  gaging  stations  (Pomeroy  upstream  and 
Huntington  downstream)  from  the  mouth  of  the  Kanawha  River. 
Pom.eroy  is  II4.  miles  above  the  tributary  while  Huntington  is 
[j.6  miles  below  it.   One  observes  that  there  is  an  upstream 
effect  of  about  1  foot  with  a  somevjhat  larger  effect  down- 
stream.  It  might  be  noted  that  the  conventional  flood  routing 
procedure  vjould  in  principle  •  furnish  an  effect  downstream 
only,  and  no  effect  upstream;  this  is  a  point  which  has  been 
discussed  at  some  length  in  the  two  earlier  reports. 
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§5.   Results  for  the  ig'-f?  Flood  at  the  Junction  of 
the  Ohio  and  Mississippi  F^ivers 

In  the  case  of  bhe  junction  problem  the  data  furnished 
vjere  better  for  our  purposes  than  the  data  for  the  Tjp-.-)er  Ohio 
River.   As  v;as  noted  above  the  data  for  tho  present  case  were 
furnished  in  such  a  v;ay  thst  average  cross  section  areas  over 
about  twenty-mile  stretches  vjere  knovm  or  easily  calculated 
and  the  resistance  coefficient  was  given  directly  over  reaches 
somewhat  shorter  than  those  in  the  Upper  Ohio,   In  fact,  only 
slight  adjustments  of  the  area  and  resistance  coefficients  had 
to  be  made;   what  these  were  will  be  described  lihortly. 

The  problera  to  be  solved  involves  stretches  each  about 
i^O  miles  in  length  along  all  three  branches  from  Cairo.   Pig. 
5.1  indicates  the  situation  schemetically.   As  boundary  con- 
ditions at  iietropolis  in  the  Ohio  River  and  at  Thebes  in  the 
Upper  Mississippi  River  the  river  stages  were  taken  from  the 
actual  records  of  the  19^4-7  flood.  (Vife  performed  another  cal- 
culation-see second  half  of  this  sectlon---in  v.'hich  the  dis- 
charge data  at  Thebes  and  Iietropolis  v/ere  used  instead  of  the 
stage  data.  )   At  Hickman  in  the  Mississippi  rlvei-  below  Cairo 
the  river  stages  of  course  were  available,  but  it  was  thouf,ht 
more  reasonable  to  make  use  of  an  average  rating  curve  at 
Hickman  as  a  boundary  condition,   V'hat  this  means  is  that  the 
effect  of  the  remainder  of  the  Mississippi  River  below  Hickman 
is  replaced  by  the  average  relation  between  discharge  and 
stage  at  Hickman  as  obtained  fx'om  flood  records.   (Figure  5.2 
displays  the  average  rating  curve  at  Hickman  and  a  few  of 
the  observed  points.)   In  addition  to  these  boundary  condit- 
ions it  is  necessary  to  if'ipose  appropriate  conditions  at 
Cairo,  the  junction.   In  effect  the  differential  equations 
are  solved  for  each  of  the  three  branches  separately  in  the 
same  manner  as  ^^/as  e::<plained  in  the  preceding  section  ror   the 
Upper  Ohio  River;  however  it  is  necessary  to  piece  together 
the  three  separate  solutions  in  the  branches  at  Cairo  by  mak- 
ing use  of  the  appropriate  continuity  conditions.   Reasonable 
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conditions  at  the  junction  are  that  the  stage  at  Cairo  should 
be  the  same  in  all  tliree  branches  and  that  the  inflow  from 
the  upper  branches  should  equal  the  flow  Into  the  Lov/er 
Mississippi.   These  conditions  together  with  the  river  stages 
and  discharges  at  all  of  the  net  points  at  the  time   t  =  0 
(v7hlch  cori'esponds  to  January  lip,  19i;7)  determine  the  flow 
uniquely  in  all  three  branches.   The  results  of  predictions 
over  a  16  day  period  are  shown  in  Pig.  5.3.   These  graphs 
furnish  the  sta<,es  as  functions  of  time  at  Cairo  and  at 
Pilckman.   As  one  sees  even  the  minor  irregularities  in  the 
observed  stages  are  followed  closely  by  the  predictions  up 
to  the  crest  of  the  flood.   The  error  vias    in  fact  never  great- 
er than  Oc6  of  a  foot.   The  amount  of   computing  time  required 
on  the  UNIVAC  for  the  16  day  prediction  was  less  than  three 
hours. 

It  is  worthwhile  to  repeat  that  the  stage  at  Hickxaan 
in  the  Lower  Mississippi  River  was  obtained  as  a  consequence 
of  assuraing  that  the  part  of  the  Jiisslsslppl  River  below 
Hickman  could  be  replaced  by  an  average  rating  curve,  that  is, 
by  an  averare  discharge-stage  relation,  at  Hickman,   The 
accurate  results  obtained  for  the  problem  in  this  way  Indicate 
that  a  large  river  system  could  be  treated  in  this  manner. 
In  other  words  even  if  oia chines  with  a  smaller  storage  capacity 
than  the  UIJIVAC  were  to  be  used  it  might  still  be  possible 
to  deal  with  lengthy  stretches  of  a  given  river  or  river  sys- 
tem by  breaking  it  up  Into  pieces  of  sufficiently  small 
lengths.   The  total  amount  of  machine  time  needed  in  comp- 
utation is  proportional  to  the  length  of  the  stretch  and  is 
not  materially  increased  by  the  above  described  subdivision. 

It  was  thought  advisable  to  carry  out  the  solution  of  the 
junction  problem  for  the  same  flood  of  19i|.7,  but  to  make  use 
of  the  known  discharges  into  the  upper  ends  of  the  Ohio  and 
Mississippi  as  a  means  of  fixing  boundary  conditions,  rather 
than  to  use  the  observed  stages  at  these  points:   In  practice. 
It  would  be  the  discharge  which  would  be  knovm  oi"*  estiinated 
rather  than  the  stage.  At  Hickman, in  the  lower  Mississippi 
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the  same  average  rating  curve,  shown  in  Fig.  5«2,  was  used 
as  for  the  earlier  calculation.   It  was  found  that  slight  ad- 
justments in  area  and  resistance  coefficients  were  needed  near 
the  upper  ends  of  the  tvjo  branches  above  Cairo  in  order  to  re- 
produce the  observed  stages  at  these  points  (i.e.  at   Thebes 
and  Metropolis)  v;hen  the  discharges  are  prescribed  rather  than 
the  stages.   The  results  with  the  new  coefficients  are  very- 
accurate,  as  one  sees  from  Pig.  S»k   which  presents  graphs 
shwoing  stages  at  Metropolis,  Thebes,  Cairo,  and  Hicknan,   Even 
the  minor  variations  in  stage  are  reproduced  faithfully,   (It 
is  seen  that  vjith  the  coefficients  used  in  the  latest  calcul- 
ation, the  stages  at  Thebes  and  Metropolis  are  reproduced 
faitlifully.   Therefore  we  would  expect  equally  good  results 
if  the  stage  upstream  were  used  as  the  boundary  data.) 

It  has  been  stated  that  proper  continuity  conditions  must 
be  imposed  at  Cairo,  the  junction  of  the  rivers,  in  order  to 
formulate  the  flow  problem  completely,   VJe  chose,  rather  nat- 
urally, the  conditions  that  the  stage  should  be  the  same  in 
all  three  branches,  and  that  inflows  and  outflows  should  bal- 
ance at  the  junction.   This,  however,  results  in  a  violation 
of  the  law  of  conservation  of  momentum,  which  is  not  very 
large,  but  does  exist.   It  might  be  more  reasonable  to  impose 
the  laws  of  conservation  of  mas.s  and  momentum,  and  put  up  with 
a  slight  discontinuity  in  stage,"  For  example,  one  could  im- 
pose the  t^^?o  conservation  lav/s,  and  require  in  addition  that 
the  stages  should  be  the  same  in  the  tv;o  upper  branches,  X'jlth 
a  slight  discontinuity  in  stage  in  the  lower  Mississippi  as  a 
result.   For  example,  in  the  model  of  the  junction  problem 
which  was  treated  in  Report  II,  a  discontinuity  of  0,2  feet  in 
stage  vjould.  result  by  making  the  calculation  in  the  manner  sug- 
gested here. 


There  would  also  be  a  slight  violation  of  the  law  of  con- 
servation of  energy,  in  that  a  loss  in  energy  would  occur. 
Such  a  loss  of  energy  at  the  junction  might  be  regarded  as 
reasonable • 
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s6.   Calculations  for  the  1950  and  19I4.8  Floods  in  Kentucky 
Reservoir 

Figure  6,1  is  a  diagrarrraatic  sketch  of  Kentucky  Reservoir 
which  shows  its  division  into  reaches  and  names  the  stations 
vjhere  the  calculated  and  observed  stages  are  compared. 

It  has  already  been  explained  in  sec.  3  l^ow  the  basic 
data  were  converted  into  the  form  needed  for  our  method  of 
computation.  We  repeat  here  that  the  data  were  accurate  for 
cross  section  areas,  since  storage  volumes  for  10-mile 
stretches  were  made  available;  and  the  resistance  coefficient 
was  essentially  ^iven  directly  by  the  basic  data.   That  the 
basic  data  were  good  is  evidenced  by  the  accuracy  with  which 
the  two  floods  were  reproduced  without  the  necessity  for  any 
but  quite  minor  changes  in  the  coefficients  once  they  had  been 
fixed  from  the  basic  data. 

In  sec,  3  it  was  already  stated  that  the  method  of 
approximating  deriv^atives  by  difference  quotients  that  x-jas 
successfully  applied  in  the  upper  Ohio  and  in  the  junction 
problems  could  not  be  used  for  the  flow  problems  in  Kentucky 
Reservoir.  The  primary  reason  for  this  was  that  the 
coefficients  A  and  B  varied  much  too  vridely  over  10-mile 
intervals,  ir^e.  at  intervals  equal  to  the  mesh  vridth.  An- 
other contributory  factor  was  the  very  rapid  increase  in  the 
discharges  into  the  reservoir  from  Pickwick  Dam.   However, 
the  difficulty  xjas  overcome  by  revising  the  scheme  of  finite 
differences  in  the  manner  explained  in  sec.  3?  while  still 
retaining  the  mesh  x^idth  of  10  miles.  We  proceed  to  discuss 
our  results  for  the  floods  of  19^0  and  19l!-8. 

In  both  floods  the  data  prescribed  at  the  tvro  boundaries, 
i,e.  at  Pickwick  Dam  upstream  and  at  Kentucky  Dam  dox^Tistream, 
consisted  in  the  observed  discharge  Q  as  a  function  of  time. 
Figure  6,2  gives  the  discharges  for  the  1950  flood.   As  one 
sees,  the  discharge  into  the  reservoir  at  Pickwick  Dam 
increases  very  rapidly  at  the  beginning,  going  from  55>000  c.f.s. 
on  Jan.  I4.  to  220,000  c.f.s.  betx^^een  Jan  6  and  7,  after  xr;hich 
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the  discharge  remains  nearly  stationary  for  about  a  week, 
after  vjhlch  there  is  a  fairly  sharp  decrease  to  about  110,000 
c.f.s.,  followed  again  by  a  slovj  rise  to  about  160,000  c.f.s. 
The  releases  of  water  from  the  reservoir  at  Kentucky  Dam 
follow  much  the  same  pattern,  but  the  changes  are  more  abrupt 
and  the  maxima  and  minima  are  farther  apart.   The  results  of 
the  numerical  solution  of  the  flov;  problem  In  the  reservoir 
are  given  in  Pig.  6,3,  which  shov;s  the  stages  at  various 
points  along  the  reservoir  as  a  function  of  time.   It  should 
be  remarked  that  the  stages  at  the  tx'jo  dams  at  the  ends  of 
the  reservoir,  as  well  as  those  at  intermediate  points,  v;ere 
obtained  by  calculation  from  the  differexice  equations.   As 
one  sees,  the  observed  flood  is  reproduced  for  21  days  very 
accurately  at  most  of  the  stations.   In  fact,  with  the 
exception  of  Savannah  and  Clifton,  where  differences  of  1  to 
1,5  feet  occurred,  the  eri'ors  are  of  the  order  of  inches. 
Upon  comparing  the  releases  at  the  dams,  as  given  by  Pig.  6.2, 
with  the  resulting  stages  given  by  Pig.  6,3,  we  observe  that 
the  large  releases  into  the  reservoir  resulted,  as  they 
should,  in  rapid  Increases  in  stage  at  the  upper  end  of  the 
reservoir.  At  Perryvllle ,  at  about  the  midpoint  of  the 
reservoir,  the  stage  inci'eased  only  slightly,  while  it 
decreased  at  Kentucky  Dam,  the  lower  end  of  the  reservoir. 
This  resulted  because  of  the  high  releases  of  water  at 
Kentucky  Dam.  The  relatively  smooth  curve  giving  the  stage 
at  Kentucky  Dam  should  be  compared  with  the  rather  wildly 
fluctuating  ci-irve  for  the  discharge  at  the  same  point.   In 
fact,  the  discharges  (and  hence  also  the  velocities)  through 
most  of  the  reservoir  in  general  vary  much  more  rapidly  than 
the  stages.  Another  comi'aent  is  also  of  interest  In  thl'? 
connection.   One  sees  that  very  heavy  releases  at  Kentucky 
Dam  seem  not  to  affect  the  accuracy  of  the  stage  calculations 
at  this  point;  the  reason  is  that  such  releases  cause  a 
receding  wave  which  lowers  the  stage  and  tends  to  smooth  out 
as  it  propagates.   On  tho  other  hand,  rapid  inflows  at 
Pickwick  do  affect  the  accuracy  of  the  approximate  soQutions, 
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since  rising  stages  are  accompanied  by  steepening  of  wave 
profiles,  and  the  finite  difference  approximations  based  on 
fixed  mesh  widths  become  less  accurate.   This  may  very  well 
be  the  cause  of  the  inaccuracies  at  Savannah  and  Clifton 
tovmrd  the  end  of  the  flow  period:   the  effect  of  the  long- 
continued  high  rate  of  discharge  at  Pickwick  Dam  has  perhaps 
led  to  inaccuracies  which  might  be  corrected  by  using  a  finer 
mesh  X'Tldth  or  by  going  over  a  more  refined  basic  scheme  of 
calculation. 

In  Pigs.  6,[|.  and  6.5  the  results  of  calculations  for 
the  I9I18  flood  are  shovm.   The  calculated  results  were 
obtained  using  the  same  coefficients  as  were  used  to  calculate 
the  flow  In  the  1950  flood.   The  results  are  given  for  7  l/k- 
days  only,  stopping  at  6:00  A.M.  on  Feb.  I3  since  the  flood 
stages  increased  from  then  on  at  such  a  rate   that  our 
approximate  method  of  computation  of  finite  differences  was 
very  Inaccurate  -  in  fact,  it  seemed  even  to  diverge.   A  finer 
mesh  width  ^x  (and  of  course,  a  shorter  time  interval  /\t ) , 
or  else  a  radical  revision  of  the  basic  method  vrould  be 
necessary  In  order  to  obtain  accurate  results  In  this  case. 
However,  for  the  first  7  l/k-   days  the  check  with  the  observed 
flood  is,  on  the  whole,  quite  good.   At  Kentucky  Dam  it  is 
very  good,  and  even  at  Pickwick  Dam,  where  the  observed  and 
computed  stages  differed  most,  the  error  is  not  excessive. 
At  all  stations,  the  shape  of  the  curve  of  the  observed  stages 
is  accurately  reproduced.   The  same  general  remarks  made 
above  for  the  1950  flood  hold  also  for  the  I9IJ.8  flood,  which 
was  also  characterized  by  heavy  releases  at  the  dams  closing 
the  two  ends  of  the  reservoir. 


The  discharge  increased  at  Plckx\'ick  Dam  from  Q  =  185,00^ 
c.f.s.  to  3)4.5,000  c.f.s.  within  l3  hours  on  Feb.  12.  Upon 
comparing  Figs.  6,3  and  6,5  one  observes  that  by  Feb.  12  the 
stages  in  the  19'j-8  flood  had  already  reached  the  maxima 
recorded  in  the  1950  flood  -  hence  fvirther  calculations  for 
the  I9I4.5  flood  could  not  lead  to  a  verification  of  the 
coefficients  based  on  the  results  for  the  1950  flood,  v;hlch 
was  another  reason  for  stopping  at  this  point. 
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The  results  of  the  calculations  for  Kentucky  Reservoir 
show,  to  our  minds,  that  the  methods  used  by  us  can  be  very 
successfullj''  used  to  solve  practical  problems  connected  v/ith 
the  operation  of  large  reservoirs.   For  example,  it  is  now 
possible  to  make  quickly  and  cheaply  a  series  of  calculations 
of  flows  through  Kentucky  Reservoir  when  various  hypotheses 
are  made  concerning  initial  states  in  the  reservoir  and  the 
discharge  schedules  at  the  txjo  ends  of  the  dam.   By  varying 
parameters  in  a  systematic  way  such  calculations  would  serve 
to  furnish  curves  useful  for  operational  purposes. 
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§7.   Discussion  of  Factors  Affecting  Accuracy^   Sug- 
gestions for  Irrnrovement  of  the  Nuraorical  Methods, 

In  carrying  out  the  coraputatlons  described  above,  insuf- 
ficient time  and  rt.oney  v;ere  available  to  improve  the  accuracy 
of  the  results  in  a  number  of  cases  in  which  there  is  little 
doubt  of  the  possibility  of  improving  them  -  often  by  rather 
obvious  devices.   Our  aim  uas  to  use  our  resources  of  time 
and  money  in  as  efficient  a  manner  as  possible  in  order  to 
carry  out  the  basic  task  of  investigating  the  general  pos- 
sibility and  practicability  of  numerical  methods  for  attacking 
flox-j  problems  in  large  rivers  and  reservoirs.   As  a  conseq- 
uence, finished  results  are  not  available  in  all  cases,  above 
all  in  the  Ohio  River,   It  v/ould  require  a  major  effort  to 
revise  the  methods  used  for  the  Ohio  River"-  in  fact,  once 
having  coded  it  in  the  wav  fii--st  thounht  reasonable,  not  too 
radical  changes  were  possible  for  us  thereafter  since  ve   wish- 
ed to  have  an  opportunity  to  study  the  other  tvjo  problems  also. 
It  seems  therefore  reasonable  to  devote  this  section  of  our 
report,  even  at  the  expense  of  occasional  repetition,  to  a 
discussion  a)  of  v;hat  our  experience  has  taught  us  about  the 
sources  of  error,  and  b)  to  offer  a  few  possibilities  for  im- 
proving the  niunerical  methods, 
a)  Factors  affecting  accuracy: 

1,   Meshwidth  /.\x,   in  general,  a  mesh  width  2/\x   =  10  miles 
seems  reasonable  in  the  three  cases  treated  by  us,  except 
in  Kentucky  Reservoir,  v.'here  the  variations  in  the  coeffic- 
ients were  too  rapid.   For  the  1950  flood  in  the  Kentucky 
Reservoir  it  was  possible  to  use  a  mesh  x^ridth  of  10  miles, 
but  only  upon  introducing  a  more  complicated  method  of 
approximating  derivatives  than  was  necessary  in  the  other 
two  cases.   For  the  19i|8  flood  in  Kentucky  Reservoir,  ten 
miles  proved  to  be  too  large  a  value  for  2/\yi   once  the  dis- 
c harges  s tarted  increasing  vjith  great  rapidity  at  Pickwick 

This  was  the  problem  with  which  vie   started.   In  retrospect, 
we  feel  now  that  it  would  have  been  better  to  have  undertaken 
it  last. 


Dam,    Probably  /\x  should  be  decreased  by  a  factor  of  tv;o 
or  three,  or,  perhaps  better,  a  different  approximation 
scheme  (perhaps  the  implicit  scheme  to  be  described  below) 
should  be  employed.   It  should  be  said,  however,  that  the 
computation  time  required  if  the  mesh  width  were  to  be  re- 
duced by  a  factor  of  two  or  tliree  would  still  be  tolerable 
in  Kentucky  Reservoir  even  if  the  Univac  were  used,  and  if 
a  faster  machine  such  as  the  IBM  70i|.  were  to  be  used,  the 
calculating,  time  would  be  negligible  even  for  these  small 
mesh  x-jldths, 

^■*      Ungated  local  inflow.   In  the  Ohio  River  the  ungaged 
local  inflow  is  a  very  Important  factor  in  the  data,  but  is 
much  less  important  in  the  other  two  problems  treated  by  us. 
In  the  Ohio  River  this  contribution  can  be  quite  high  -  in 
the  reach  St.  Hary«  s  -  Pomeroy,  for  example,  it  is  k.0  %   of 
the  total  flow,  and  in  many  cases  it  is  as  much  as  one  and 
a  half  times  the  gaged  inflow  from  tributaries.   In  the 
xrhole  stretch  from  V'heeling  to  Cincinnati,  there  is  one  time 
when  the  discharge  of  the  main  stream  at  V/heeling  is 
160,000  cfs  while  /^.yo^OOO  cfs  flows|over  the  banks.   The 
accuracy  with  which  the  ungaged  flow  has  been  given  is,  of 
course,  not  precisely  known,   We  suspect  it  to  be  not  very 
accurate.   For  example,  a  recent  check"'"'  of  the  inflow  from 
the  drainage  areas  shows  that  two  different  methods  of  using 
the  unit  hydrograph  method  for  converting  rainfall  into  dis- 
charge (the  difference  consisted  in  taking  6-hour  rather 
than  12-hour  time  intervals)  resulted  in  discharges  which 
differed  by  five  to  ten  thousand  cfs  for  a  period  of  half  a 
day. 

The  authors  of  this  report  feel  that  the  inaccuracies 
in  calculating  the  19'+8  flood  in  the  Ohio  may  well  stem  in 
large  part  from  inaccuracies  in  the  local  inflows.   The  final 
results  obtained  for  the  19i;5  flood  were  good,  but  they  were 


This  may  also  have  been  a  factor  in  affecting  the  accuracy  in 
theJ)hio  River,  where  the  m.aximum  rate  of  Increase  of  sta4  was 
"t  ~  ^'^  ftVhr»,  \^hich  is  even  larger  than  the  maximum  of  0.75 
ft,/hr,  in  Kentucky  Reservoir. 

Made  while  this  report  v/as  being  written. 
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brought  about  by  somewhat  radical  changes  in  cross  section 
area  and  resistance  coeff icients.   If  it  is  true  that  such 
changes  i-jere  made  to  compensate  for  erx'-ors  in  inflov/  rather 
than  for  inaccuracies  in  coefficients,  it  would  net  be  sur- 
prising that  subsequent  calc.ilr.t ions  for  another  flood  based 
on  such  coefficients  vjould  prove  to  be  Inaccurate.   At  the 
same  time,  it  should  also  be  pointed  out  that  there  are 
still  other  3D  urces  of  ei'ror  in  the  Ohio  River  problem,  which 
were  not  investigated  -  above  all  the  effect  of  a  mesh 
VJidth  of  10  miles  in  the  upper  sections  where  the  arsa  co- 
efficients varied  considerably,  and  vjhere  rates  of  increase 
in  sta^e  were  larger  than  those  vh  ich  occurred  in  Kentucky 
Reservoir. 
b)   Suggestions  for  improving  the  accuracy: 

1  •   Different  interpolation  schemes  for  cross  s ect^ion  ai'eas. 

In  sec,  3  our  method  of  dealinc  with  cross  section 
areas  was  explained  in  details   For  the  purpose  of  the  dis- 
cussion here  the  main  po3.nt  ic   that  average  cross  sections 
were  obtained  for  each  of  the  reaches  (in  the  Ohio  these 
vjere  60  -  90  miles  long),  these  values  were  taken  to  hold 
at  the  centers  of  the  reaches,  and  values  at  intermediate 
net  points  were  obtained  by  linear  interpolation  from  the 
center  of  one  reach  to  the  center  of  the  next  adjacent  reach. 
The  resulting:  cross  section  area  at  a  ^iven  place  then  de- 
viates often  quite  widely  from  the  local  value,  especially 
at  the  gaging  stations  themselves  which  are  usually  placed 
at  narrow  parts  of  a  river.   However,  our  results  show  that 
the  stages  obtained  by  calculations  based  on  such  average 
cross  section  areas  are  accurate  locally,  but  that  the  dis- 
charges obtained  by  multiplying  the  velocities  obtained  by 
us  into  the  cross  section  areas  used  in  our  calculations  can 
be  quite  wrong  locally.   This  point  has  already  been  discuss- 
ed in  section  k   in  connection  x-^/ith  the  Ohio  River  problem, 
and  a   satisfactory   means  to  obtain  correct  discharges 
was  given.   Nevertheless,  it  v/ould  probably  be  better  to 
revise  our  scheme  of  interpolation  for  cross  section  areas 


56 

and  resistance  coefficients  in  such  a  way  as  to  bring  the 
model  of  the  river  v/hich  the  averaging  process  yields  into 
somexjhat  closer  agreement  vrith  the  actual  river.   One  vjay  to 
do  so  would  be  to  take  the  actual  cross  section  areas  at  the 
ends  of  the  reaches,  then  employ  the  total  storage  volunes 
for  the  reaches  in  such  a  vjay  as  to  obtain  an  average  cross 
section  at  the  center  of  the  reach,  so  that  linear  interpol- 
ation between  it  and  the  actual  areas  at  the  ends  of  the 
reach  would  yield  the  correct  total  storage  volume  in  the 
reach.   One  could  expect  that  correct  values  for  the  discharge 
at  the  gaging  stations  would  then  be  obtained  simply  by  mul- 
tiplying the  calculated  velocity  with  the  cross  section  area, 
which  nox^f  vjould  be  the  correct  local  value,   V/e  carried  out 
this  program  partially  for  the  Ohio  River,  but  obtained  bad 
results  because  x-.'e  did  not  recalculate  resistance  coeffic- 
ients, but  used  the  same  values  as  before*   It  is  hovrever, 
clear  that  the  local  value  of  G   should  be  U5ed  at  the  gag- 
ing stations,  since  the  average  value  used  in  the  earlier 
calculation  can  depart  widely  from  the  local  value,  as  we 
know  from  sec.  I4.,      Time  did  not  permit  a  revision  of  G   in 
an  appropriate  fashion,  but  vie  vjould  expect  an  improvement 
in  the  results  if  it  were  done.   In  addition,  this  new  viay 
of  dealing  va  th  area  and  resistance  coefficients  vjould  not  be 
essentially  harder  to  apply  nor  lead  to  more  complications 
in  coding  than  the  former  scheme. 
2 e   Smoothing  of  cross  s ecj: j.qn  data. 

In  Kentucky  Reservoii'  (cf,  sec,  3  and  the  preceding 
sec,  6)  we  have  seen  that  a  new  scheme  for  approximating  der- 
ivatives was  made  necessary  because  of  the  rapid  variation 
in  the  area  coefficient  along  the  reservoir.   The  data  for 
cross  section  areas  v;ere  supplied  as  averages  for  reaches  ten 
miles  in  length  (obtained  from  topographic  maps),  and  these 
values  differed  by  factors  as  large  as  two  in  adjacent 
reaches.   Nevertheless,  the  data  in  this  form  \-jere  used  by 
us,  but  a  revision  in  the  numerical  methods  was  found  to  be 
necessary.   It  might  have  been  better  to  have  smoothed  the 
average  cross  sections  first,  and  then  employed  the  smoothed 
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values  as  coefficients  in  the  differencial  equations.   A- 
gain,  time  did  not  permit  a  trial  of  such  a  method, 
2 0   Cross  section  are£. s  from  topopraphic  ni&ps  for  the  Upoer 
Ohio  Rivero 
In  view  of  the  good  results  obtslned  in  the  junction 
problem  and  in  ICentucky  Reservoir,  in  vrhich  avera>-,e  cross 
section  areas  were  obtained  from  topographic  maps,  it 
xTOuld  probably  be  advisable  to  do  the  same  thing  for  the 
Upper  Ohio  River.   There,  it  vjill  be  recalled,  the  average 
cross  sections  v;ere  obtained  from  storage  volumes  in  long 
reaches  and  thevSe  in  turn  resulted  from  balancing  out  in- 
flows and  outflo^^rs  in  actual  floods.   Thus  the  cross-section 
areas  used  by  us  contained  errors  due  to  faulty  inflow  data- 
and,  as  we  have  seen,  the  ungaged  inflow,  which  can  only 
be  estimated  with  not  too  good  accuracy,  can  be  a  quite 
considerable  fraction  of  the  total  inflow,   ive  would  re- 
comiaend  in  general  the  use  of  average  cross  sections  obtain- 
ed from  maps  since  the  extra  labor  involved  would  not  be  a 
large  fraction  of  the  total,  and  one   vjould  thus  be  sure  of 
the  accuracy  of  one  of  the  most  essential  basic  quantities, 
^'   Implicit  schemes  for  numerical  solution  of  the  differen- 
tial equations. 
It  has  been  stated  repeatedly  that  our  method  of  num- 
erical solution  of  the  differential  equations,  which  pro- 
ceeds by  advancing-  the  values  of  the  stage  and  discharge 
from,  knovm  values  at  time  t  to  new  and  as  yet  unknown 
values  at  the  time  t  +  /^  t ,  requires  that  the  time  step 
/^^t  must  be  kept  very  small  (of  the  order  of  9  minutes  in 
our  thiree  cases)  in  comparison  with  the  total  time  for  which 
flow  calculations  ivould  normally  be  made  (of  the  order  of 
weeks).   The  time  required  for  making  computations  on  a  com- 
puting machine  is  inversely  proportional  to  the  time  inter- 
val /\t.   In  our  case,  it  did  not  matter  too  much,  since 
the  amount  of  time  needed  for  calculation  on  the  Unlvac  v;as 
quite  reasonable.   Nevertheless,  if  much  lonuer  stretches 
of  a  river  or  river  system  were  to  be  treated  numerically. 
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this  factor  of  machine  tine  mi.iht  become  crucial.   For  this 
reason  (and  for  others  to  be  mentioned  later),  it  would  be 
worthwhile  to  study  a  radically  different  method  of  numerical 
solution  in  whicli  the  values  of  stage  and  velocity  can  be  ad- 
vanced in  time  steps  of  the  order  of  hours  instead  of  minutes, 
though  at  the  expense  of  a  more  complicated  com.putational 
scheme. 

VJe  proceed  to  discuss  briefly  one  possible  way  to  set   up 
a  numerical  scheme  of  solution  which  might  permit  values  for 
/,t   of  t\'o  or  three  hours.   In  the  accompanying  figure  7.1 
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a  double  rovj  of  net  points  in  the  x,t-plane  is  shown,  and 
the  object  of  the  calculation  is  to  advance  the  values  of  v 
and  H  usinj:  all  of  the  values  in  both  rows  simultaneously. 
To  this  end  the  differential  equations  (2.1)  with  A,  replac- 
ing BH.  are  written  down  for  each  i-ect angle,  e.g.  LRUP: 

(7.1)  — 2/:^ ^  - — \ — sT^x 


%,vl  +   Gpv2  +  G,v2  +  Gpv| 


'U^U 


'pvp 
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(Hp  +  Hp)  -  (Hy  +  H^) 


Qv 


U,?,L,R 


(7.2) -^^3, ^~ 

The  unknowns  In  these  equations  are  v„,  Vp  H^,  and  Hp,  I.e. 

two  more  than  there  are  equations.   However,  if  one  boundary 
condition  is  given  at  each  of  the  upstream  and  do^^;nstreara 
ends  (discharge,  say,  or  a  ratin<;^  curve),  the  number  of 
equations  and  the  numbor  of  unknox-.rns  would  be  the  same.   In 
Pig,  7.1,  for  example,  there  v.-ould  be  eight  values  to  be 
fixed  at  the  four  upper  net  points;   there  ;-J0uld  be  six 
equations  of  the  type  (7«1)  and  (7.2)  and  these  together  x^?lth 
two  values  at  the  boundaries  (or  two  relations  there),  vjould 
yield  a  system  of  equations  the  same  in  number  as  the  number 
of  unknovms,"  However,  in  marked  contrast  to  the  situation 
in  the  schemes  used  by  us  so  far,  the  present  system  of 
equations  would  be  a  system  of  nonlinear  simultaneous  equat- 
ions which  could  be  solvable  practical.ly  only  by  methods  of 
successive  approximation.   Probably  the  best  method  to  use 
for  solving  the  equations  would  be  an  iterative  method,  with 
Initial  trial  values  for  the  unknox-jns  selected  by  linear  ex- 
trapolation fror-i  the  results  at  previous  rows  of  net  points. 
Since  we  know  that  v  and  H  do  not  chan^'e  much  for  tliaes 
as  great  as  2  or  3  hours,  it  seems  quite  possible  that  such 
a  method  would  yield  good  approximations  to  the  exact  solut- 
ion v;hen  time  steps  of  these  mat:nitud.es  are  used. 


"It  is  the  implicit  character  of  this  scheme  that  differen- 
tiates it  from  the  previous  schemes.   The  fact  that  all  of  the 
values  at  all  net  points  for  a  time  t  are  used  s irault anoou ily 
in  determining  the  values  for   t  +  /\t   is  the  reason  why  the 

slope  of  the  characteristics  no  longer  dictates  the  maximum 
size  permissible  for  /\t. 
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The  implicit  scheme  would  have  still  another  advantage 
over  former  schemes  which  mirht  Increase  the  accuracy  of  the 
approximation.   The  advantage  is  that  the  irnplicit  sciieme  pro- 
vides automatically  a  means  to  ensure  that  the  continuity  con- 
dition holds  in  the  large,  i.e.  that  Inflows  and  outflows 
vjould  balance  not  only  at  each  net  point,  but  over-all,  if 
such  a  coraputation  vjere  to  be  made.   The  methods  used  bv  us 
hitherto  do  not  always  check  well  for  a  whole  river  -  in  fact, 
in  the  Ohio  fairly  large  discrepancies  are  found  upon  ruaking 
a  check  for  the  whole  river  (This  is  understandable  siiice 
the  computed  discharges  have  to  be  adjusted  as  described  at  the 
end  of  section  q.).   The  sane  remark  would  also  apply  to  the 
lav;  of  conservation  of  momentum,  thou'h  with  less  force  since 
we  find  this  law  to  be  well  satisfied  in  the  large  when  our 
methods  are  used. 

Such  implicit  schemes  have  been  proposed  before  for  sim- 
ilar problems,  but  as  far  as  v/e  know,  they  have  never  been 
used  on  a  large  scale.   It  is  also  of  interest  to  observe  that 
H.  A,  Thomas  (The  Hydraulics  of  Flood  Movements  in  Rivers, 
Carnegie  Institute  of  Technology, , 1937 )  proposed  a  variety  of 
schemes  for  numerical  solution  of  river  wave  problems,  some  of 
which  would  probably  diverge,  but  also  amon/'  them  is  a  scheme 
essentially  the  same  as  the  imr)liclt  scheme  proposed  here. 

^»   Simplified  differential  equatiqn,'j_^ 

It  might  be  possible  to  simplify  the  calculations  some- 
what by  simplifying  the  differential  equations,  either  through 
the  omission  of  terras  felt  to  be  unimportant,  oi'  by  develop- 
ments and  expansions  of  various  kinds  wtiich  yield  simpler 
equations. 

For  example,  the  term  vv   in  the  dynamical  equation  could 
probably  be  omitted  without  causing  much  error,  and  also  the 
term  qv/A  (the  contribution  of  the  inflow  to  the  momentum). 
However,  as  long  as  one  operates  numerically  with  a  system,  of 
hyperbolic  equations  (or  even  parabolic  equations)  the  savings 
from  simplifications  in  the  differential  equations  are  not 
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very  great  since  the  great  bulk  of  the  v;ork  consists  in  the 
processing  of  the  basic  data,  and  coding  it.   Still,  It  is 
by  no  means  impossible  that  further  studies  could  lead  to 
formulations  in  which  less  elaborate  calculating  equipment 
than  the  Unlvac  would  suffice,  at  least  for  some  types  of 
problems.   For  example,  a  flood  in  a  long  river  for  wiiich  the 
inflovjs  occur  mainly  in  tiie  head-vjaters  so  that  the  flood  pro- 
pagates essentially  only  in  the  downstream  direction,  v/ith 
little  or  no  backwater  effects,  very  likely  could  be  treated 
by  simpler  methods  than  those  used  by  us  in  this  report. 
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§8.   Model  Studies  Compared  with  Numerical  Studies  Using 
Digital  Computers 

In  the  preceding  sections  some  renarks  have  already  been 
made  concerning  the  relation  between  the  methods  explained 
here  for  the  computation  of  flood  waves  in  rivers  through  use 
of  the  basic  differential  equations  as  compared  with  the  method 
of  using  physical  models  of  a  river  or  river  system.   These 
remarks  will  be  amplified  somewhat  in  this  section. 

A  brief  description  of  the  procedures  used  in  m.aking 
model  studies  should  first  be  given.   Such  models  -  at  least 
those  models  of  the  Mississippi  system  which  are  novj  at  the 
Waterways  Exr^eriment  Station  at  Jackson,  Mississippi  -  are 
built  on  a  rather  large  scale,  covering  acres  of  ground  in 
fact.  They  are  built  vjith  a  very  much  exaggerated  vertical 
scale  compared  with  the  horizontal  scale.   Models  are  made  of 
concrete  vrith  the  channel  being  built  up  accurately  from 
topographical  surveys;  this  is  a  costl3'-  and  time  consuming 
feature.  The  inflows  must  be  fed  into  the  channel  by  rather 
complicated  machines  (pumps)  which  can  reproduce  any  given 
discharges  as  functions  of  time.   However,  because  of  cost 
these  are  not  placed  at  all  of  the  main  tributaries.   Instead, 
the  inflows  over  a  considerable  stretch  of  the  river  are 
lumped  together  to  form  a  composite  hydrograph;  this  hydrograph 
is  reproduced  by  a  machine,  and  the  resulting  discharge  is  fed 
into  the  river  at  a  number  of  points  in  various  proportions. 
The  actual  inflows  are  thus  not  put  in  accurately  at  the 
points  X'jhere  they  enter  the  river,  but  rather,  a  certain 
average  inflow  is  put  in.   (It  is,  however,  probably  a 
reasonable  average . ) 

The  river  stages  at  the  main  gaging  stations  are  then 
recorded  electronically  and  the  records  can  be  seen  in  a 
special  house  containing  recorders  for  each  of  the  gaging 
stations.  Air  conditioning  is  necessary  (just  as  vjith  the 
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UNIVAC  computer)  because  of  the  large  anount  of  electronic 
equipment  that  is  involved.   All  of  this  equipment  is  quite 
costly,  and  its  operation  reqviires  a  considerable  staff. 

It  has  already  been  mentioned  that  one  of  the  main 
physical  featiires  governing  the  flow  in  a  river,  namely  the 
roughness  of  the  river  bed  which  gives  the  flow  its  turbulent 
character,  cannot  be  scaled  properly  for  a  model.   Instead 
it  is  necessary  to  introduce  obstacles  in  the  bed  of  the  model 
in  order  to  reproduce  the  floods  actually  observed.  At 
Jackson,  these  obstacles  take  the  form  of  small  brass  knobs 
screwed  into  the  bed  of  the  model,  or  of  wire  screen,  the 
latter  placed  as  a  rule  in  shallov/er  parts  of  the  stream. 
The  fact  is  that  the  resistance  force  is,  together  with 
gravity,  the  biggest  force  conditioning  the  flov;  in  the  stream 
and  the  fact  that  this  element  of  the  problem  cannot  be  scaled 
is  a  good  reason  for  calling  such  models  analogue  computing 
machines  rather  than  true  scale  models. 

It  has  been  stated  above  that  the  process  of  num^erical 
computation  of  floods  is  analogous  in  all  respects  to  the 
method  of  predicting  floods  using  a  model.  The  only  difference 
is  that  we  have  made  use  of  average  cross  section  areas  and 
resistance  coefficients  rather  than  actual  local  values  as  is 
done  in  the  models  .   It  is  the  requirement  of  reproducing 
highly  accurate  cross  sections  that  constitutes  one  of  the 
costly  features  in  building  models.  The  experience  reported 
above  for  the  case  of  the  junction  problem  and  for  Kentucky 
Reservoir  shows  that  it  is  not  necessary  to  use  actual  cross 
sections  but  rather  that  it  suffices  to  use  average  cross 
sections  over  10-20  mile  stretches.  Even  in  the  Ohio  River, 
where  averages  for  reaches  of  60  to  9^  miles  in  length  were 


On  the  basis  of  our  experience  using  average  cross  sections 
it  is  even  indicated  that  models  might  also  be  built  using 
average  cross  sections  ratlier  than  actual  cross  sections,  with 
the  expectation  of  getting  results  having  sufficient  accuracy, 
and  presumably  at  smaller  cost.  This  in  fact  is  an  idea  which 
via.3   advocated  long  ago  by  the  vjell-knovm  hydraulics  engineer 
H.  A.  Thomas. 
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used,  the  fact  that  the  results  were  not  as  pood  as  in  the 
other  two  cases  was  quite  likely  less  due  to  the  use  of 
averages  over  such  lone  stretches  than  it  was  to  the  fact  that 
the  basic  data  did  not  furnish  sufficiently  accurate  averages. 

The  next  vital  element  involved  both  in  the  model  and 
in  the  method  of  numerical  computation  concerns  correct 
determination  of  the  resistance  factor.   In  models  this  is 
done  in  a  completely  empirical  way  in  the  manner  indicated 
above.   In  making  our  numerical  studies  vje  could  at  least 
begin  with  a  first  reasonable  estimate  for  this  quantity  based 
on  analysis  of  actual  observations.   And,  in  the  case  of  the 
jxinction  problem  and  Kentucky  Reservoir,  where  the  area  and 
resistance  coefficients  were  well-determined  by  good  basic 
data,  little  or  no  further  experimentation  with  the  resistance 
coefficient  was  required  to  reproduce  the  flood.  True,  it  ims 
found  necessary  in  the  less  accurately  defined  Ohio  River 
problem  to  correct  the  resistance  factor  by  making  predictions 
for  a  relatively  short  period  and  comparing  with  an  actual 
flood.   This  corresponds  exactly  to  what  must  be  done  in  a 
model  vjhen  brass  knobs  or  wj.re  screen  are  added  or  taken  ai\ray 
from  various  parts  of  the  model  in  order  to  reproduce  observed 
floods. 

There  are  various  ways  in  which  the  method  of  numerical 
computation  is  much  more  flexible  than  the  method  of  using 
models.   In  the  first  place  it  is  quite  easy  to  vary  cross 
section  areas  if  it  is  thought  desirable  to  improve  the 
accuracy.   In  a  model  that  cannot  be  done  readily;  it  perhaps 
should  be  done  in  some  Instances  especially  vjhen  the 
topographical  data  which  were  used  for  constructing  the  model 
are  very  old  and  there  are  possibilities  that  the  river 
channel  had  changed  in  the  meantime.  Another  feature  in  which 
numerical  methods  are  more  flexible  concerns  the  method  of 
dealing  with  the  local  runoff  and  the  flow  from  tributaries. 
In  operating  with  models  it  vras  indicated  above  that  the 
inflows  are  fed  into  the  model  at  a  relatively  small  number 
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of  points,  after  composite  unit  hydrographs  have  been 
calculated,  since  the  apparatus  needed  for  this  purpose  is 
rather  complicated.   In  making  numerical  studies  very  little 
additional  effort  is  required  to  introduce  the  flows  at  any 
point  that  might  be  desired.   Of  course  the  method  of  finite 
differences  requires  tliat  the  fl.ows  be  Intr'oduced  in  the 
intervals  between  the  net  points.  This,  however,  means  in 
the  cases  studied  by  us  that  flows  are  introduced  at  points 
only  10  miles  apart.  That  is,  in  each  10  mile  interval  tne 
flow  either  from  a  tributary  entering  that  interval  or  from 
the  local  runoff  is  Introduced. 

It  is  sometimes  argued  that  the  method  of  using  a  model 
has  an  advantage  over  the  method  of  numerical  calculation 
since  vrith  the  model  it  is  possible  to  observe  what  happens, 
say,  in  case  a  dike  breaks  and  a  flow  out  of  the  river  channel 
takes  place.   However,  it  would  seem  to  us  that  the  results 
of  such  observations  would  be  Illusory  unless  observations  on 
floods  of  this  kind  were  available  for  the  purpose  o£   making 
verification  runs  in  the  model;  vjhlle  on  the  other  hand,  if 
such  data  were  available  they  could  be  used  to  carry  out 
predictions  by  numerical  methods.   In  fact,  as  we  have  seen, 
it  would  be  possible  to  estimate  reasonable  values  for  the 
resistance  coefficient  T^jithout  any  flow  data,  and  proceed  hj 
numerical  calculation  to  study  flows  where  no  such  data  were 
available. 

In  the  end,  presumably,  it  is  the  m.atter  of  cost  which 
is  the  dominant  feature  in  any  comparison  of  the  two  methods. 
So  far  onljr  the  limited  experience  of  this  one  group  Is 
available  for  estimating  costs  of  the  numerical  method. 
However,  in  three  years  a  relatively  small  group  has  been  able 
to  carry  out  the  work  outlines  above  starting  with  no  previous 
experience.  The  total  amount  spent  for  the  purpose  of 
research  and  development  of  numerical  methods  including  the 
cost  of  using  the  WJIVAC  was  well  under  $200,000,00.   With  the 
experience  gained  now  it  should  be  possible  to  code  the  data 
for  equivalent  stretches  of  other  rivers  for  much  less  than 
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this.  As  for  the  amount  of  machine  time  needed  for  actual 
computation,  we  have  seen  it  to  be  quite  small.   For  the 
junction  problem,  for  example,  a  l6  day  prediction  required 
less  than  3  hours.  The  rental  fee  for  the  UNIVAC  from  the 
Army  Map  Service  in  I/ashington  is  only  $75.00  an  hour; 
privately  oxmed  calculating  equipment  is  usually  rented  at 
somewhat  higher  rates  than  this  but  even  so  this  cost  is  not 
large.   if  a  model  is  built,  it  and  the  equipment  which  goes 
with  it  must  be  serviced  and  maintained,  one  would  think  at 
rather  high  cost.   Once  the  data  for  a  river  have  been  coded 
for  a  machine,  it  is  stored  permanently  in  a  few  rolls  of 
magnetic  tape  and  is  ready  at  any  moment  to  be  used  to  solve 
any  flow  problem;  codes  would  also  be  prepared  in  advance  so 
that  the  special  initial  and  boundary  data,  and  the  inflows, 
all  of  ^^Jhich  differ  from  case  to  case,  could  be  rapidly 
prepared  for  use  in  the  machine. 

Finally,  it  should  be  remarked  once  more  that 
calculi  ting  equipment  is  constantly  being  improved,  so  that 
our  experience  based  on  the  use  of  a  UNIVAC  is  not  at  all 
final.   In  fact,  the  use  of  the  IB?4  TOli.  would  probably  reduce 
the  calculating  time  needed  for  the  problems  di:-cussed  in 
this  report  in  the  ratio  of  about  1  to  1^.   Thus  even  the 
calculating  time  for  the  most  difficult  case,  the  Ohio,  vjould 
be  reduced  from  hours  (for  flows  lasting  vreeks)  to  a  matter 
of  minutes. 

Thus  the  authors  are  of  the  or)inion  that  numerical 
methods  offer  great  advantages  in  comparison  with  model 
studies  for  the  types  of  oroblems  dealt  vjith  in  this  report 
-  i.e.  wave  problems  in  long  rivers  and  reservoirs.   However, 
this  should  not  be  taken  to  imr)ly  any  adverse  criticism  of 
model  studies  in  general,  not  even  of  those  made  in  the  past 
for  large  rivers  -  after  all,  the  kind  of  computing  eouipment, 
and  knowledge  of  numerical  analysis,  that  was  used  for  the 
problems  discussed  in  this  report  has  been  available  only  for 
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a  few  years.   Only  in  cases  like  the  present  ones,  in  uhich 
a  reliable  mathematical  formulation  of  the  problems  in  a  not 
too  complicated  Xvay  is  possible,  can  model  studies  be 
dispensed  with.   On  the  other  hand,  model  studies  could 
hardly  be  dispensed  \-jith   for  such  a  problem  as,  for  example, 
the  determination  of  the  flox^r  characteristics  through  tlie 
gates  and  tributaries  of  Kentucky  Dam  itself,  since  an 
adequate  purely  mathematical  formulation  of  such  a  problem 
would  be,  if  not  dovniright  impossible,  then  at  least 
enormously  complicated. 
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19.   Conversion  of  Rainfall  Data  into  Rionoff  Data  for  the 
Ohio  River 

The  progress  of  a  flood  wave  may  depend  decisively  on 
the  amount  of  water  which  is  contributed  by  the  ungaged  runoff 
of  rainwater  coming  from  the  area  drained  by  the  river:   that 
is,  the  flow  from  the  main  valley  and  from  small  tributaries 
which  are  not  gaged.   The  simple  arithmetical  calculation 
which  converts  rainfall  data  into  runoff  data  has  been  coded 
and  tested  on  the  UNIVAC.   It  took  less  than  10  minutes  of 
machine  time  to  compute  the  runoff  data  for  the  5  reaches  for 
a  10  day  period.  All  of  our  flood  predictions  viere  made  with 
runoff  data  that  was  hand-computed  and  supplied  to  us  by  the 
Corps  of  Engineers.  VJe  merely  report  here  on  the  obvious  fact 
that  it  is  feasible  to  prepare  the  runoff  data  on  the  UWIVAC 
and  indicate  the  method  that  was  used. 

The  Ohio  River  Division  of  the  Corps  of  Engineers 
supplied  us  with  the  data  on  rainfall  over  each  of  the  5 
drainage  areas  which  are  associated  with  the  5  reaches  into 
which  the  upper  Ohio  River  is  divided.   That  is,  the  rainfall 
amounts  for  each  12  hour  period,  from  S:00  P.M.  on  Feb.  25, 
I9I4.5  to  midnight  of  Mar.  6,  were  given  together  with  the 
percentages  of  these  which  represent  the  amount  of  water 
(called  excess  rainfall)  which  flows  into  the  river  according 
to  certain  fixed  unit  hydrograph  proportions.   Our  calculations 
were  made  for  the  period  from  Feb,  26,  1911-5  through  Mar.  6, 
19i|5«   (There  ^^?as  no  rain  for  i;  days  prior  to  Feb,  26.) 
VJe  used  equation  (9.1)  to  convert  the  12  hour  rainfall  data 
for  a  reach  into  runoff  data  for  the  reach  at  6  hour  intervals 
(since  the  basic  interval  into  which  we  divide  the  UNIVAC 
computation  is  6  hours  -  that  is,  we  represented  q(x,t)  as 
being  a  linear  function  of  the  time,  t,  for  a  period  of  6  hours). 
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ru2.r(t-36)+u.  .r(t-2l+)+u^.r(t-l2)+UQ'r(t)  , 


for  t  =  6:00  A.H.  or 
6:00  P. 14. 


(9.1)   q(t)  = 


I  u^.r(t-[|.2)+u--r(t-30)+u^.r(t-l8) 


+u  .r(t-6)+u  .r(t+6)  , 


I 


for  t  =  noon  or 
midnight 


where  q(t)  is  the  amount  of  overbank  inflow  for  the  reach 
expressed  in  units  of  one  thousand  cubic  feet  per  second, 
vrhile  r(y)  is  the  araount  of  excess  rainfal]  in  inches  that  fell 
over  the  drainage  area  during  the  12  hour  period  preceding  y, 
i.e.  from  (y-12)  until  y;  and  u-.  jUp, . . .  ,Ug  (see  Table  9.1)  are 
the  6  hour  instantaneous  hydrograph  coefficients  for  the 
drainage  area.   Note  that  only  the  rainfall  of  the  'oreceding 
2  1/2  days  is  permitted  to  contribute  to  the  runoff 


"Beach 

VJheeling 
St.   Marys 

St.    Marys 
Pomeroy 

Pomeroy 
Huntington 

Huntington 
Haysville 

Haysville 
Cincinnati 

^^1 

11,172 

311,557 

30,727 

Ll.2,328 

13,969 

^2 

1)4, '+57 

1-14,721 

39,761; 

51+,  778 

18,078 

^3 

l6,J+29 

50,820 

1+5,187 

62,21+7 

20,51;3 

\ 

18,^00 

56,918 

50,609 

69,717 

23,008 

u^ 

20,536 

63,525 

56,1+83 

77,309 

25,679 

^6 

17,71+3 

51;,  885 

1+6,801 

67,227 

22,186 

^7 

12,650 

39,131 

31].,  79'+ 

1;7,930 

15,818 

^3 

10,51'-! 

32,525 

28,919 

39,838 

13,ll;7 

j^ 

7,229 

22,361 

19,882 

27,389 

9 ,  039 
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The  inflow  data  that  v;as  used  in  our  Ohio  River  flood 
calculations  was  read  off  at  6  hour  intervals  from  curves  of 
inflow  prepared  by  the  engineers  x^^ith  the  use  of  12  hour 
hydrograph  coefficients.  He  compare  the  values  of  inflow 
obtained  by  us  through  the  use  of  6  hour  hydrograph 
coefficients  in  equation  (9.1)  with  the  data  supplied  by  the 
engineers;.  Pig.  9.1  is  for  the  reach  Wheeling  to  St.  Marys 
and  Pig.  9.2  is  for  the  reach  St.  Marys  to  Pomeroy.   It  is 
seen  that,  as  is  to  be  expected,  the  6  hour  coefficients 
produced  slightly  higher  crests  and  lower  troughs  than  was  the 
case  for  the  12  hour  coefficients  (since  the  peaks  didn't  occur 
at  the  ends  of  the  12  hour  periods). 

In  particular,  we  noted  that  the  short  dry  period  xi/hlch 
occurred  on  March  5>  coincides  with  a  brief  drop  in  stage  in 
the  measured  hydrograph  at  the  stations  between  Wheeling  and 
Cincinnati  on  that  day.   Our  calculated  hydrographs  do  not 
have  as  pronounced  a  dip  and  v;e  attribute  this  fact  in  p^.rt 
to  our  having  used  the  12  hourly  smoothed  runoff  data  which 
do  not  have  as  pronounced  a  dry  period  on  March  5  as  actually 
occurred. 
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COMPARISON     OF    RUNOFF    INFLOWS     OBTAINED    WITH 
6  HOUR   AND   12    HOUR    COEFFICIENTS 

(WHEELING— ST.  MARYS) 
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Fig.  9. 1 
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COMPARISON    OF    RUNOFF    INFLOWS    OBTAINED 
WITH    6    HOUR    AND     12    HOUR     COEFFICIENTS 

(ST.  MARYS -POMEROY) 
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Fig.  9.2 
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